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ABSTRACT 


The Numerical Electromagnetics Code (NEC) was used to 
evaluate the admittance, average power gain, and the electric 
near and far field of a monopole antenna mounted on a cubical 
box over a perfectly conducting ground plane. Two models of 


the box, employing surface patches and wire grids, were 


evaluated. The monopole was positioned at the center, the 
edge, and at a corner of the box's top surface. Admittance 
and average power gain of the antenna were calculated. NEC 


results were examined and compared with experimental data and 


with results frome! eaten. another independent 
electromagnetic modeling code. The near electric field was 
calculated for both models. Computer graphics techniques 


were presented for plotting NEC near field results using 
DISSPLA (Display Integrated Software System and Plotting 
LAnguage), a commercial graphics package. Contour and 328 
amplitude, and phase plots of the near electric fields were 
presented. Radiation patterns were calculated to relate far 
field and near field behavior of the antenna. Surface patch 


and wire grid models are compared and conclusions were 


presented. 
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LSS INTRODUCTION 


A. GENERAL BACKGROUND 

The arrangement of HF antennas aboard ships requires the 
consideration of many factors. Among these are: input 
impedance, gain, blockage angles, patterns, weight, firing 
arcs of weapons, and transmission line length restrictions. 
Additional factors that greatly influence naval operations 
are RF radiation hazards (RADHAZ) and electromagnetic 
interference (EMI). 

Shipboard operations are carried out within fixed (small) 
distances from HF transmitting antennas, so the Navy has a 
unique and long-standing operational problem--the radiation 
from these antennas can be hazardous to personnel, ordnance, 
fuel, and electronic equipment because of the intensity of 
the fields close to the radiating elements (near fielday. 
The goal of the antenna designer is to select a location for 
the antenna that permits it to perform with a minimum of 
restrictions placed on personnel activity, electronic 
communication equipment, fuels and ordnance. 

In the selection of a site for an antenna, it becomes 
evident that all of the desired requirements cannot be met, 
therefore, trade-offs have to be made based on a priority 
system determined by the ship's primary mission. For 
instance, on a combatant ship, placing antennas relative to 
gun mounts, missile launchers, and magazines are Of  Magor 


Z 


importance. On an auxiliary, however, the underway 
replenishment gear might receive the greatest attention. 

For the above reasons, the Navy has recognized the need 
and focused its interest on the investigation of the near 
fields of antennas for a number of years. 

Near field structure is very complex, and previous 
theoretical analysis is practical for only very simple 
antennas in uncomplicated geometrical settings. With the 
advent of high-speed computers, approximate solution 
techniques such as the Method of Moments became practical. 
This technique is suitable for the calculation of 
electromagnetic fields anywhere, including in close proximity 
to the radiating element. The degree of accuracy of the 
solution is a function of the number of computations required 
to obtain that solution (computer-time). 

The Method of Moments technique is the theoretical basis 
for the Numerical Electromagnetics Code (NEC), which is a 
code for the simulation and analysis of the electromagnetic 
response of antennas and other metal structures [Ref. 1]. 
NEC is the computer simulation tool that was used in this 


investigation of near fields. 


B. REASONS FOR THIS STUDY 
1. Hazards of Electromagnetic Radiation (RADHAZ) 
Since near fields of Navy HF antennas are primarily 
responsible for radiation hazards (RADHAZ), the accuracy of 
near field predictions must be adequate to ensure that 


Z 





personnel and equipment can continue to function safely in 
modern shipboard environments. There are three types of 
radiation hazards: 


a. Hazards of Electromagnetic Radiation to Personnel 
(HERP) . 


b. Hazards of Electromagnetic Radiation to Ordnance 
(HERO) . 


c. Hazards of Electromagnetic Radiation to Fuels 
(HERF) . 


These hazards and their general near field criteria are 
discussed in detail in Reference 2. 
2. Electromagnetic Interference (EMT) 

Near fields of shipboard antennas can also be the 
cause of radio frequency interference. On board ships RF 
interference problems can be quite complex. [In _ land 
installations, there are several preventive measures for 
reducing interference, such as isolation of communication 
related equipment and separation of transmitting and 
receiving antennas to avoid near field interaction. Ona 
ship, it is almost impossible to do this. The design is 
constrained by the nature of the ship and care has to be 
taken in operating radiating equipment in order to minimize 
interferences. 

Because shipboard transmitting and receiving antennas 
are close together and the ship's superstructure provides a 
very large nonlinear surface to interact with communications 


Signals (near field interaction), intermodulation problems 


are likely because of currents induced by nonlinear devices 
in the near field region of an antenna. 

Intermodulation sources are nonlinear items present 
in the immediate surroundings of the receive or transmit 
antennas. Metallic life lines are major sources of 
interference; expansion joints and mooring or anchor chains 
used to secure the ship at the pier are found to introduce 
significant nonlinearities. Other nonlinear items include 
ladders, guard rails, antenna guying wires, and booms. [Ref. 
3] 

These environmental nonlinear sources can be 
classified as two major categories, (1) nonlinear resieamge 
"Junctions" created by metallic mating surfaces, and (2) 
ferromagnetic materials. Ferromagnetic materials, such as 
steel, generate interfering frequencies when driven by strong 
fundamental currents caused by the near fields of an antenna. 

Table 1 from Reference 4 is a listing of a number of 
Significant losses in lives, fleet assets, and dollars due to 


RADHAZ and due to Electromagnetic Interference (EMI). 





TABLE 1 [Ref. 4] 
MAJOR LOSSES DUE TO EMI OR EMI "FIXES" 


During Time of Military Active Environment: 
-- 32 aircraft/134 crew lost, $172M damage, EMI triggered A/C 
rocket detonation on USS Forrestal flight deck - Late 60's. 


-- Most of the topside destroyed, crew injuries, $100+M 
damage, anti-radiation missile detected friendly cruiser - 
Late 60's. 


-- All equipment lost, crew fatalities, $10+M damage, anti- 
radiation missile detected USMC radar - Early 70's. 


-- Frigate lost, crew fatalities, $200+M replacement, HMS 
Sheffield secured ASMD radar to use SATCOM to avoid EMI-1984. 


During Peacetime Environment - Late 70's thru 1986: 

-- Planes Crashed: NATO aircraft (A/C) overseas - A/C lost, 
crew fatalities, $10+M replacement. USN aircraft in 
California - A/C lost, crew fatalities, $10+M replacement. 


-- FFG ASW helo while hovering over FFG, ship radiation 
caused blade servo fault, helo damage, crew injuries, $750K 
damage. 


Ship Missiles Lost and Near Miss Hits of USN Ships: 
-- Cruiser destroyed missile to avoid another cruiser, 1 
missile and 1 target drone lost, $10+M. 


Ship Damage Due to EMI: 
-- Terrier launcher rail fire, target drone lost, $100K+ 
replacement. 


Seecrart dropped from foilborne operation due to ship control 
system EMI - Bow damaged, crew injuries unknown, $10+M. 


-- Fires in propulsion room due to EMI "fixes" installed to 
prevent total shutdown of plant in presence of RF, minor 
damage, no injuries, crew repaired. 


Equipment Damage Due to EM Environment: 
-- ECM equipment burned out by radiator on neighbor ship, 
equipment lost, $1M+ replacement. 


Sailors RF Burned: 
Common occurrence for crew. Increases noted from 1982-1985: 


-- Personnel RADHAZ - ten-fold increase 
-- Personnel RF burns - six-fold increase 
-- Hazards to ordnance - five-fold increase 


Cc. SCOPE OF THE sluUDy 

The Numerical Electromagnetics Code (NEC) is used to 
evaluate the effects of the shipboard topside environment on 
the admittance and the electric-near field structure of a 
6 cm monopole antenna mounted on a cubical box. 

A ship's superstructure has been approximated by two 
theoretical scaled models of cube-shaped boxes of 0.1 m sides 
(A/3 at a frequency of 1 GHZ) over a perfectly conducting 
ground plane. For the first model, the surface patch 
modeling technique is used and for the second, the wire grid 
modeling technique. 

The antenna is A/5 at a frequency of 1 GHZ and is mounted 
at the center, edge, and corner of the top surface. NEC 
results for admittance are obtained for these three 
positions. 

For the surface patch model, the results are compared to 
measured experimental data from References 5 and 6. For the 
same model, NEC results are compared with "PATCH" code 
results. "PATCH" is another independent code recently 
developed using a much different theoretical and numerical 
algorithm implementation [Ref. 7]. The computational 
accuracy of NEC is then examined and discussed. For the wire 
grid model, the admittance characteristic results from 
Reference 8 are compared with the experimental data mentioned 
above. Comparison between surface patch and wire grid 


numerical modeling techniques are made. 


NEC is then used for calculation of near electric fields 
for surface patch and wire grid models. Fertran programs 
were constructed for plotting NEC results of near electric 
fields in conjunction with DISSPLA (Display Integrated 
Software System and Plotting LAnguage), [Ref. 9], which is a 
commercial proprietary plotting package. 

Contour plots and 3-D plots of near electric field are 
presented. The phase variation of the near electric field is 
examined and plotted also. The effect of the shipboard 
topside environment on the electric near field of the 
monopole antenna is discussed. 

For the surface patch model radiation patterns were 
calculated in order to relate the far field and near field 


behavior of the antenna. 


D. THEORETICAL BACKGROUND 
1. Near Field Theory 
The fields around an antenna may be divided into two 
regions, one near the antenna called the near field or 
Fresnel zone and one at a large distance called the far- 


field or Fraunhofer zone. Figure la indicates these two 


regions. 
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Figure la. Fields of a Radiating Antenna [Ref. 8] 

The boundary between the far field and near field is 
the distance r = 2D°/4 , where D is the maximum length of the 
antenna in meters and A is the wavelength in meters. The 
distance from the surface of the antenna to this boundary is 
called the near field region, while beyond this boundary the 
region is called far field. The near field region is divided 
into two subregions -- the reactive and radiating near field. 
The reactive near field is extended to A/2n from the 
antenna's surface, while in practice a distance of A 
represents this boundary. The phase of the magnetic and 
electric field is almost quadrature in regions within a 
wavelength of the antenna (reactive near field). Beyond the 


distance of a wavelength, the electric and magnetic fields 





are propagating in phase (radiating near field) until the far 
field is reached. In the far field, the shape of the field 
pattern is independent of the distance, while in the near 
field the shape depends on it. 

Consider a dipole antenna which is enclosed in an 
imaginary sphere. Figure 1b shows the energy flow near the 


dipole. 
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Figure 1b. Energy Flow Near a Dipole Antenna [Ref. 8] 


The waves are reflected near the poles while in the equator 
are diffusing perpendicular to the dipole. Consequently, 
energy is oscillating near the antenna while part of it flows 
outward in the region of the equator. The oscillating energy 
represents reactive power, while the outward flow represents 
radiated power. The near field is related to charge and 
current density while the far field is related only to 


Surtrent . 


oe Numerical Electromagnetics Code (NEC) - Method of 


Moments 

NEC is a computer code which analyzes the 
electromagnetic reponse of antennas. The program uses an 
electric field integral equation (EFIE) and a magnetic field 
integral equation (MFIE) to do this analysis. The EFIE and 
MFIE integral equations are solved in NEC, for current 
distribution using the method of moments approach. This 
approach reduces the integral equation to a system of linear 
algebraic equations in the form of a matrix, where the 
unknowns are coefficients in an expansion of the current. 
The matrix equations can be solved for the current in a high- 
speed computer. There are two limitations to the moment 
method: (1) the amount of computer storage necessary for the 
elements of the impedance matrix and (2) the time required to 
calculate those elements and solve the system of equations. 
When the current distribution is computed, all effects of the 
antenna structure can be determined. [Ref. 1] 

During this investigation, two modeling approaches 
are used, surface patch modeling and wire grid modeling. In 
surface patch modeling a two-dimensional field equation and 
appropriate surface boundary conditions are applied to solve 
for two orthogonal components of current on small patches 
representing the actual surface. Either an electric field 
integral equation (EFIE) or a magnetic field integral 
equation (MFIE) may be used theoretically. NEC employs the 
MFIE method while the code mentioned previously, "PATCH", 
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uses an EFIE method. The wire grid modeling approach 
approximates the antennas and the surfaces on the ship with 
thin wires using the EFIE method. Surfaces are represented 
as wire grids in which the grid cells are much smaller than 
the wavelength. Either modeling approach can provide data 
upon which engineering judgments can be _ made. The 
investigations of the two modeling approaches in the past 
shows that wire gridding provides slightly more accurate data 
in comparison to surface patch modeling. This is another 
point of investigation of this thesis. However, surface 
patch modeling can apparently offer some amount of savings in 


computer time cost. [Ref. 2] 


ee 


IIT. SURFACE PATCH MODELING 


A. GUIDELINES FOR MODELING USING NEC 

A conducting surface is modeled by means of multiple, 
small flat surface patches, corresponding to segments used to 
model wires. The patches are chosen to completely cover the 
surface to be modeled. The parameters defining the surface 
patch are the cartesian coordinates of the patch center, the 
components of the outward-directed, unit normal vector and 
the patch area. The program computes the surface current on 
each patch along the orthogonal unit vectors which are 
Cangent to the surface. Four -patch shape options exist: 
arbitrary shape, rectangular, triangular, and quadrilateral. 
For the rectangular, triangular, and quadrilateral patehece 
the outward normal vector ff is specified by ordering of 
corners 1, 2, and 3 using the right-hand rule. The vectors 


A Aw 
t, are parallel to the side from corner 1 to corner 2, andues 


NN 
=x t,. The coordinates of corners 1, 2, and 3 have to be 
specified in the progran. Figure 2 shows how the multiple 
patches are defined for a rectangular surface. Rectangular 


patches were used for the modeling purposes of this thesis. 


eZ 


NY 





Figure 2. Rectangular Surface Covered by Multiple Patches 
[Ref. 1] 

When a wire is connected to a surface, the wire must end 
at the center of a patch with identical coordinates used for 
the wire end and the patch center. The program then divides 
the patch into four equal patches about the wire end. Only a 
Single wire may connect to a given patch and a segment may 
have a patch connection on only one of its ends. A wire may 
never connect to a patch formed by subdividing another patch 
for a previous connection. The use of surface patches is 
restricted to modeling voluminous bodies. Very long narrow 
patches should be avoided when subdividing the surface. [Ref. 


1] 


B. LIMITATIONS AND OPTIMIZATION CRITERIA 

The accuracy of the results depends on the patch size 
measured in wavelengths. A minimum of about 25 patches 
should be used per square wavelength of surface area, with 
the maximum size for an individual patch less than or equal 
to 0.04 square wavelengths. 


as) 


In NPGNEC (an NPS customized version of NEC-3) which was 
used for simulation purposes in this thesis, the number of 
patches plus the number of segments must be less than 300. 
When the total number is greater than 300 and less than 
1,000, the NPG1000 version of NEC-3 must be used; this is 
done for the wire grid modeling in the second part of this 
thesis. Since the program does not integrate over patches 
except at a wire connection (active patch), the patch shape 
does not affect the results. The parameters affecting the 
results are the location of the patch centroid, the patch 
area, and the outward unit normal vector. 

The optimum model for complex structures which the first 
part of the thesis investigated can be estimated by varying 
the segment and patch density, observing the results and the 
convergence of the solution. In the case of an edge- 
mounted antenna, the accuracy of the results are expected to 
depend upon the size of the patches near the edge. Smaller 
patches are suggested at edge areas since the current 


magnitude may vary rapidly in this region. 


C. PREVIOUS STUDIES 

In References 5 and 6, experimental and computational 
investigations were performed to determine the admittance 
characteristics of a monopole antenna mounted on a cubical 
conducting box of 0.1 m sides (473 at frequency 1 GHz) overa 
ground plane. The antenna, a 6 cm monopole (4/5 at frequency 
1 GHZ), was tested for three different mounting positions on 
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the top surface. Experimental admittance data for the 6 cm 

monopole antenna were plotted vs. frequency [Ref. 5]. 

Numerically calculated data using the "PATCH" computer code 

were also tabulated [Ref. 5]. "PATCH" is a frequency domain 

electromagnetic analysis code based on a method of moments 
solution to the Electric Field Integral Equation (EFIE) [Ref. 

Piles In this code objects are modeled by planar triangular 

patches which easily conform to surfaces and boundaries of 

general shape and allow variable patch densities over the 
surface of the object. 

Di. MODEL OF 6 CM MONOPOLE ANTENNA MOUNTED ON A CUBICAL 
CONDUCTING BOX--ADMITTANCE AND AVERAGE GAIN INVESTIGATION 
Monopole antennas mounted on arbitrary shaped objects 

have received much attention because they represent practical 
situations where antennas are mounted on aircraft, ships, 
vehicles, and buildings. The antennas are not always 
centrally mounted and are often located near edges and 
corners. 

Super-structures of Navy ships can be modeled as 
conducting boxes which form a "ship-like" structure. Navy 
shipboard topside antennas are often monopoles and are 
attached to the top of these boxes, in different locations. 

The numerical model of a cubical five-sided box of 

0.1 m per side was constructed using NEC (the bottom was not 

included as a surface since the box was placed on a perfectly 


conducting ground plane). A 6 cm monopole antenna was placed 


1S 


at the center, at the edge (3.63 cm from Center) anawaa, 
corner (5.14 cm on a diagonal from center), as shown in 
Figure 3. This model configuration was selected since 
experimental data was available for it for comparison from 


Reference 5. 
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Figure 3. Geometry of Monopole Attached to the 
Conducting Box Over a Ground Plane 

The first part of the investigation checked the input 
impedance and average gain of the model using NEC-3 (NPGNEC). 
Comparisons of calculated values for input impedance and 
average gain were made to experimental values from Reference 
5 in order to validate the computer model before proceeding 
with an investigation of near fields. 

This modeling exercise compared driving point admittance 


from NEC as patch density is varied, to both experimental 
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values and PATCH code results from Reference 5. The number 
of patches on the top of the box was varied in search of an 
optimum value of surface samples, which would later be used 
for near field calculations. The top was subdivided to 
retain symmetry, as much as possible, and to closely match 
positions of the antenna on the experimental model. The 
monopole was positioned at the center, at an edge (3.63 cm 
from the center) and at a corner (5.14 cm on a diagonal from 
the center), very close to the experimental model positions. 

Nine frequencies from 1 to 1.4 GHZ were selected. NEC 
input data files for the three different locations of the 
monopole are presented in Appendix A. 

The Numerical Green's Function option (NGF) was used to 
exploit partial symmetry, reducing solution time. The 
nonsymmetric part, the 6 cm monopole antenna, was added to 
the NGF file of the box for impedance and average gain 
Calculations. 

1. Monopole at the Center 

The monopole was divided into five segments and 
placed at center of the top surface. The top surface was 
Sivided into 25 (5x5), 49 (7x7), 81 (9x9), and 121 (11x11) 
patches. The impedance (in OHMS) and admittance (in mMHOS) 
are listed in Tables 4, 5, 6, and 7 for the nine frequencies. 
Experimental data from Reference 5 are shown in Table 2. NEC 
statistics for the four top surface configurations are in 


Table 3. The number of patches used for each case, the patch 


Is 


area, and the variation of run time are shown. Double 
precision NEC (DNPGNEC) was used for the 5x5 and 7x7 patch 
case to investigate the possibility of better accuracy than 
single precision NEC. The results are in Tables 8 and 9. 

To gain insight into the accuracy of the results, 
the average power gain was calculated. Correct average power 
gain does not necessarily insure accurate results. 
Therefore, input impedance was also checked. The average 
power gain provides a good check of the solution accuracy if 
it is close to two for radiating systems located over 


perfectly conducting ground which is our case, or close to 


one for free space conditions. Average power gain within 
10% of theoretical is adequate for engineering purposes. For 


all cases, the calculated average power gain is between 1.96 
(minimum) to 2.02 (maximum) which is in good agreement with 
theory. 

NEC results and measurements, for the monopole at the 
center, while the patch density varies on the top surface 
(Case 1: 5x5, Case 2: 7x7, Case 3: 9x9, and Case 4: lixiije 
are presented in Figures 4, 6, 8, and 10. In Figures 5, 7, 
9, and 11, the PATCH code results have been superimposed on 
those of NEC and the measurements for comparison. The PATCH 
code is compared to experimental results in Figure 12. 

Single precision (NPGNEC) calculations vs double 
precision (DNPGNEC) results (Tables 4, 6, 8, and 9) for 


conductance and susceptance vary only in the second and 
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third decimal places. This difference is negligible and 
cannot be observed ina plot. For all cases, the antenna was 
fed at the base segment. 

CASE 1: 5x5 patches on top surface (Figures 4 and 5) 

The correlation between measurements and NEC 
calculations for admittance in Figure 4, was generally quite 
poor. Apparently, the solution has not converged. There is 
good correlation for susceptance between 1.2 and 1.4 GHz. 
PATCH results of Figure 5 correlate well with measurements in 
comparison to NEC. 

CASE 2: 7X7 patches on top surface (Figures 6 and 7) 

NEC values of admittance correlate well with 
measurements (Figure 6), but the run time observed in Table 3 
is almost double compared to Case 1. The solution has 
apparently converged. The increased accuracy is probably due 
to the top of the surface patch box having been more finely 
divided into smaller patches which provides finer details of 
surface currents. Observing Figure 7, we see that NEC has 
somewhat better correlation with measurements than PATCH 
does. 

CASE 3: 9x9 patches on top surface (Figures 8 and 9) 

There is good correlation of NEC results with 
measurements in Figure 8, especially in conductance, but not 
quite as good as in Case 2. 

NEC performs better than PATCH for both conductance 


and susceptance. As the number of patches increases (81 


be 


patches on the top surface vs 49 for Case 2 and 25 for Case 
1), the run time is doubled and tripled correspondingly. 
CASE 4: 11x11 patches on the top surface (Figures 10 and 11) 
The correlation of NEC results is good especially in 
conductance (Figure 10) and is better than PATCH (Figure 11). 
Run time is doubled compared to Case 3, as can be observed 
from=Table 32 
GENERAL RESULTS: The solution seems to converge at Case 2 
with 7x7 = 49 patches on the top surface, patch area 0.00020 
m2 (0.0022 square wavelengths at 1 GHz) and 0.00005 m? 
(0.00055 square wavelengths at 1 GHz) for each of the four 
central "subpatches" where NEC subdivides a normal size patch 
at the connection point to the monopole. NEC gives excellent 
results compared to measurements for Cases 2, 3, and 4 and 
has better performance than PATCH for these three cases. 
PATCH results are better than NEC results for Case 1 (5x5 = 
25 patches on the top). For easier comparison between the 
two codes, NEC vs. PATCH admittance values for the Cases 1, 
2, 3 and 4 are plotted separately in Figures 13, 14, 15, and 
16. In Figures 17 and 18, NEC's conductance and susceptance 
vS Measurements as the top surface patch density increases 
(5x5, 7x7, 9x9, and 11x11 patches) are available jiea 
comparison. Consequently, the optimum model appears to be 
Case 2 (7x7 = 49 patches on the top), the one that is 
selected for the near field investigation for surface patch 


numerical model. 
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TABLE 2 
MEASURED DATA FOR A 6CM MONOPOLE ON A CUBICAL BOX 


6CM MONOPOLE POSITION 
ADMITTANCE (mMHOS) 


FROM CENTER FROM CENTER 
. 2 0 ieee 
0 iLO. 7 
3 2 














FREQUENCY 


(GHZ) 










0 . 
aid) : 
1 7 0 : 
1 iO . 2 21.0 . 
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TABLE 3 
NEC STATISTICS FOR MONOPOLE AT CENTER OF PATCH BOX 
AS THE NUMBER OF PATCHES ON TOP INCREASES 
(NPS MAINFRAME COMPUTER - IBM SYSTEM 370) 


NUMBER OF PATCH RUN TIME (sec) 
PATCHES AREA (m2) |SINGLE PRECISION NEC 


surface on sides 
18 x 4 = 12) 02000565 
i937 
surface on top 
1 5x5 == 25} 0.00040 
4 











monopole attach- 
ment patch 0.00010 6.747 


ova. vearones soa 
surface on sides 

25-863 
surface on top 
2 7x7= 49} 0.00020 








monopole attach- 
ment patch 4| "0.00005 8.940 









0.00056 
54.170 


0.00012 


surface on sides 
18 x 4 = 72 
surface on top 
9x9= 8 
3 
monopole attach- 
ment patch 0.00003 1307 








il 

4 
surface on sides 
18 x 4 = Le 0.00056 
107.00 
surface on top 
LY xs i. = 21 0.00008 


monopole attach- 
ment patch 4} 0.00002 a7 27 67 
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TABLE 4 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(5 X 5 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
POWER 


ADMITTANCE (mMHOS) 
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TABLE 5 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(7 X 7 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
POWER 
GAIN ADMITTANCE (mMHOS) 
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TABLE 6 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(9 X 9 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
POWER 
GAIN ADMITTANCE (mMHOS) 
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TABLE 7 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(11 X 11 PATCHES ON THE TOP SURFACE = BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
(GHZ) POWER 


ADMITTANCE (mMHOS) 
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TABLE 8 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(5 X 5 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 
(DOUBLE PRECISION CALCULATION) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
(GHZ) POWER 


ADMITTANCE (mMHOS) 
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TABLE 9 
CALCULATED DATA 
MONOPOLE AT CENTER OF SURFACE PATCH BOX 
(7 X 7 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 
(DOUBLE PRECISION CALCULATION) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
(GHZ) POWER 


ADMITTANCE (mMHOS) 
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Figure 4. Monopole at Center of Patch Box 
(5x5 Patches on Top) 
NEC Admittance vs Measurements 
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Monopole at Center of Patch Box 
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NEC Admittance. vs Measurements and "PATCH" Code 
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Figure 5. 
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Monopole at Center of Patch Box 


Figure 6. 


(7x7 Patches on Top) 
NEC Admittance vs Measurements 
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Figure 9. 
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Monopole at Center of Patch Box 


(11x11 Patches on Top) 
NEC Admittance vs Measurements 


Figure 10. 
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Figure 11. Monopole at Center of Patch Box 
(11x11 Patches on Top) 
NEC Admittance vs Measurements and 
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Figure 13. Monopole at Center of Patch Box 
(5x5 Patches on Top) 
NEC Admittance vs "PATCH" Code 
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Figure 14. Monopole at Center of Patch Box 
(7x7 Patches on Top) 
NEC Admittance vs "PATCH" Code 
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Figure 16. Monopole at Center of Patch Box 
(11x11 Patches on Top) 
NEC Admittance vs "PATCH" Code 
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Figure 17. 
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NEC Susceptance vs Patch Density 


Figure 18. 
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2. Monopole at Edge 

The 6 cm monopole was divided into five segments and 
placed at an edge (Fig. 3). The measured data [Ref. 5] are 
for a monopole attached to an edge at a distance 3.5 cm from 
the center. Two computer model configurations of the top 
surface were selected in an attempt to match the geometry of 
the experimental model: 81 (9x9) patches which corresponds to 
a monopole position of 3.33 cm from the center and 121 
(11x11) patches which corresponds to 3.63 cm. The difference 
in distance for the position of the monopole in the NEC model 
compared to the actual physical geometry, is 0.17 cm less and 
0.13 cm more, correspondingly. 

For each configuration, the frequency was varied from 
1 to 1.4 GHz in nine computer runs. The impedance (in OHMS), 
admittance (in mMHOS) and calculated values of the average 
gain are shown in Tables 10 and 11. DNPGNEC was used for 
the case of 9x9 patches on top, in order to investigate the 
possibility of higher accuracy than single precision NEC; the 
results are shown in Table 12. From Tables 10, 11, and itz 
the calculated average gain varies from 1.89 (minimum) to 
1.92 (maximum), which is in good agreement with the 
theoretical value of 2.0. 

Results of NEC predictions vs. measurements for the 
monopole at the edge, with top surface patch density varying 
(Case 1: 9x9, Case 2: 11x11), are presented in Figures 19 and 


21. Figures 20 and 22, show PATCH code results of Reference 
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5 superimposed with NEC calculations and include measurement 
plots for comparison. PATCH code and experimental results 
are presented in Figure 23. Experimental data for the 
monopole located at the edge are shown in Table 2. The total 
number of patches and the patch area are listed in Table 3. 
To compare performance between NEC and PATCH for each case, 
see Figures 24 and 25. 

Comparing NPGNEC vs. DNPGNEC results (Tables 11 and 
2), we observe that the values of conductance and 
susceptance differ only in the second and third decimal 
place. This difference is negligible and cannot be observed 
in plots of NPGNEC vs DNPGNEC results. For all cases, the 
antenna was fed at the base segment. 

CASE 1: 9x9 patches on the top surface (Figures 19 and 20) 

The agreement is very good for conductance (G) and 
susceptance (B), when NEC results are correlated with 
measurements (Figure 19). 

In Figure 20, a conductance comparison between the 
performance of NEC and PATCH vs. measurements shows 
performance to be almost the same between frequencies of 1 to 
tls GHZ. 

CASE 2: 11x11 patches on the top surface (Figures 21 and 22) 

For this configuration, the position of the monopole 
in the NEC model is closer to the experimental position (3.63 
cm vs 3.5 cm) than in the 9x9 patch case. NEC conductance 


appears to follow the same general trend as measurements. 


45 


PATCH's conductance values are 17% closer to 
measurements than NEC. NEC and PATCH have almost identical 
performance for susceptance and both have the same resonant 
frequency of 1.20 GHz vs 1.16 GHz for themeasurement. 

GENERAL RESULTS: 

In both cases, NEC resultS are reasonably close to 
experimental data. For clarity, results for NEC and PATCH 
for Case 1 are plotted in Figure 24 and those for Case 2 are 
plotted in Figure 25. 

NEC performs better than PATCH in Case 1 and vice 
versa for Case 2. The solution seems to converge at Case 1 
(9x9 patehes ontop). However, the model of Case 2 is 
selected for the near field calculations because the position 


of the monopole is closer to the experimental model position. 
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TABLE 10 
CALCULATED DATA 
MONOPOLE AT EDGE OF SURFACE PATCH BOX (3.33CM FROM CENTER) 
(9 X 9 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
POWER 
GAIN ADMITTANCE (mMHOS) 
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TABLE 11 
CALCULATED DATA 
MONOPOLE AT EDGE OF SURFACE PATCH BOX (3.63CM FROM CENTER) 
(11 X 11 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
(GHZ) POWER 
GAIN ADMITTANCE (mMHOS) 
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TABLE 12 
CALCULATED DATA 
MONOPOLE AT EDGE OF SURFACE PATCH BOX (3.63CM FROM CENTER) 
Crt X 11 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 
(DOUBLE PRECISION CALCULATION) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
(GHZ) POWER 
GAIN ADMITTANCE (mMHOS) 
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(9x9 Patches on Top) 
NEC Admittance vs Measurements 


Figure 19. 
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3. Monopole at Corner 
The 6 cm monopole in the NEC model was placed at the 


corner, 5.14 cm on the diagonal from the center and fed at 


the base. The position of the monopole for the experimental 
model is 5.15 cn. The top surface was divided into 121 
(11x11) patches. The calculated values for the impedance, 


admittance and average power gain are shown in Table 13. 
Experimental data are presented in Table 2 and the total 
number of patches as well as the patch area in Table 3. 
Average power gain is between 1.80 (minimum) and 1.83 
(maximum) and is within the 10% tolerance of the theoretical 
value. 

NEC results of admittance vs. measurements are 
presented in Figure 26. PATCH code results from the 
calculated data of Reference 5 have been superimposed on 
NEC's data, for comparison, in Figure 27. PATCH code results 
are compared to experimental results Figure 28. NEC vs. 
PATCH is compared in Figure 29. 

NEC is in excellent agreement with measurements in 
both conductance and susceptance (Figure 26). PATCH seems to 
perform slightly better than NEC (Figure 27), however, 
conductance and susceptance in the range of 1.15 to 1.40 GHZ 


are virtually identical for both codes. 
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TABLE SEs 
CALCULATED DATA 
MONOPOLE AT CORNER OF SURFACE PATCH BOX 
(5.14CM ON THE DIAGONAL FROM CENTER) 
(11 X 11 PATCHES ON THE TOP SURFACE - BASE SEGMENT FEED) 


FREQUENCY AVERAGE IMPEDANCE (OHMS) 
POWER 
GAIN ADMITTANCE (mMHOS) 
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Eee NEAR ELECTRIC FIELD 

Near fields are more difficult to calculate than far 
fields. When calculating radiation in close proximity to an 
antenna, the antenna no longer appears as an infinitesimal 
point (as is the assumption for far field), and the energy 
radiating from it comes from separated locations insofar as 
the strength of radiation from each sub-section of the 
radiator is concerned. This means that terms in the field 
expressions with powers of 1/r* (r is the distance from the 
origin of the antenna to the field point) are appreciable in 
magnitude compared to the 1/r dependent terms which are 
dominant in the far field. The complex Poynting vector 
m2 (E x H*) thus contain terms with higher powers of (1/r), 
in addition to the radiation field term. The near field is 
thus very dependent on the charge density and the current 
while the far field is affected only by the current. 

For near field calculations, NPGNEC computes’ the 
magnitude of the peak electric field (E-Total) in (V/m), 
which is the vector sum of the three components Ey, Ey, and 
Ez (the calculation of Epeak iS a special feature that has 
been added to the NPGNEC version of NEC which is in use at 
NPS) . NEC outputs the magnitude and phase of each component 


Ex, Ey, and E, separately. 
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1. Near Flectric Field in 2=Dimensions for tnelolen 
Monopole Antenna on the Cubical Conducting Box 


Using the optimum models in NEC for the monopole at 


the center (7x7 patches), at the edge (11xil patches) and the 


corner (11x11 patches), we proceed with the near field 
investigation. The NE card is used in conjunction with the 
PL card; both are described in detail in Reference 1. The 


input data files are presented in Appendix A for the near 
field calculations and for the three positions of the 
monopole on the patch box. For all cases, NEC was tasked to 
produce two hundred points in x, y, or z directions, with a 
spatial variation of 0.01 m. The starting point for Ghesueau., 
field calculation is 0.05,A away from the wire and box metal 
structure. The calculation of the near field is extended to 
a distance of at least 6),A in the x, y, or z directions, 
respectively. In that manner, we ensure that the field 
points are far enough away from any metallic structure and at 
a distance where the near field variation can be observed in 
detail. The plots of the near field are presented in Figures 
31 and 32 (Monopole at center -- variation in x and 2), 
Figures 33, 34, and 35 (Monopole at edge -- variation in x, 
y, and z), and Figures 36, 37, and 38 (Monopole at corner-- 
Variation in xX; VY, andes 

In order to compare NEC near fields for the monopole 
on the box with theoretical values, we consider a linear 
current element I=IpeJwt of length Az oriented in the z 
direction and located at the origin as in Figure 30 (Ref. 
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nO) - For convenience, assume Ip is the amplitude of the 


current. This antenna is a simple radiating structure but it 


will demonstrate basic properties of the near electric field 


for all small linear antennas. 


N 








Figure 30. A Linear Current Radiator 


The complete electric-field intensity of the antenna 
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The only part of the field entering into the expresstommionr 
the radiated power is that part consisting of the terms 


varying as Cone that is; 


This is called the radiation field. The parts of the field 
varying as r~? and r~? are called the induction field. The 
induction field does not represent an outward flow of power, 
but instead depicts a storage of reactive energy in the 
vicinity of the radiating current element (Reference 10). 
Consequently, the terms that are functions of distance r of 


the E-Field in Eq. 2.1 are: 


All the other terms are phase terms or constants. Generally, 
the magnitude of electric field has an r dependence which can 


be expressed as: 


(2525 
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For the monopole at the box center, the near electric field 


is presented in Figure 31. From this figure and NEC output, 
we observe for a distance close to zero, r = 0.0166 m, the 
mmide of the field is |E(r)| = 45.9 (V/m). Using this value 


in Eq. 2.2, we can find the constant factor C = 0.000222065 
(normalization factor). Using the calculated value of the 
constant C in Eq. 2.2 and varying the distance r from the 
linear current element, we calculate the corresponding |E(r) | 
for each distance; the r-dependence is plotted in Figure 39. 
Comparing Figure 31, for the near electric field of the 
monopole on the box center with the near field of the current 
element in Figure 39, both near electric fields follow the 
same trend which is: a maximum in close proximity to the 


antenna with a gradual reduction as the observation point 


moves away. A knee of the curve appears approximately at 
0.125 m (0.41A away from the antenna). The results for the 
monopole at the edge and corner are almost the same. The 


"spikes" observed in Figures 31, 33, and 34 are from the edge 


of the box. 
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MAGNITUDE VS DISTANCE eraign 
monopole at center of patch box 7x7 patches on top 





Magnitude Ee aK (v/m) 
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Figure 31. Monopole at Center of Patch Box 


(7x7 Patches on Top) 
Near Electric Field on xX Axis 
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Magnitude EpeaK (v/m) 


MAGNITUDE VS DISTANCE freg=Ighz 
monopole at center of patch box 7x7 miches on top 
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Figure 32. Monopole at Center of Patch Box 
Near Electric Field on Z Axis 
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Magnitude EpeaK (v/m) 


MAGNITUDE VS DISTANCE freq Ighz 
near e-field 6 cm monopole at edge 3.63 cm from center 
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Figure 33. Monopole at Edge of Patch Box 
Near Electric on X Axis 
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MAGNITUDE VS DISTANCE freq | ghz 
near e-field 6 cm monopole at edge 3.63 cm from center 
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Figure 34. Monopole at Edge of Patch Box 
Near Electric Field on Y Axis 
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MAGNITUDE VS DISTANCE freq | ghz 
near e-field 6 cm monopole at edge 3.63 cm from center 
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Figure 35. Monopole at Edge of Patch Box 
Near Electric Field on Z Axis 
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Magnitude EpeaK (v/m) 








MAGNITUDE VS DISTANCE pes Ighz 
em on diagonal 


near e-field 6 cm monopole at corner 5.1 
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Figure 36. Monopole at Corner of Patch Box 
Near Electric Field on X Axis 
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Magnitude EDEAK (v/m) 


MAGNITUDE VS DISTANCE ned Ighz 
near e-field 6 cm monopole at corner 5.14 cm on diagonal 
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Figure 37. Monopole at Corner of Patch Box 
Near Electric Field on Y Axis 
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Near Electric Field of a Current Element 


Figure 39. 
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Ps NEAR ELECTRIC FIELD MAGNITUDE AND PHASE CONTOUR PLOTS 
Fortran programs were developed to plot the magnitude and 
phase of the near electric field using NEC output data. The 
programs use subroutines of the DISSPLA (Display Integrated 
Software System and Plotting LAnguage) graphics software 
package [Ref. 9]. Appendix B lists the Fortran source codes 
and describes their function. PROGRAM I: Contour plot data 
sets of the magnitude Epeak aS well as the magnitude of the 
Paetividual components E, and E, of the electric field (a 
section of a plane in 3-dimensional space), composed of a 
30 X 30 point window in space, are produced. PROGRAM ITI: 
This is similar to Program I but for the phase of x and z 
components of the electric field, for 100X100 points. 
PROGRAM III: MThree-dimensional plot data sets of Epeak are 
produced. For these plots, a 3-D map of the electric field 
is shown on the top and the 2-D contours of the field are 
presented on the bottom. CONTOR FORTRAN: This is a 


subroutine which is called by the above three programs and 


creates the plots. Positive value contours are presented by 
solid lines and negative contours by dashed lines. In these 
programs, NEC output for the near electric field is 


interfaced with the plotting facilities of DISSPLA. 

For reference, two simple antennas with known results are 
examined: (1) a 0.15 m (A/2) dipole antenna in free space, 
and (2) a 0.06 m (A/5) monopole antenna over a perfect ground 


plane. NEC input data files which produce magnitude and 


Wa 


phase of near electric field for the two antennas are listed 
in Appendix B. 


pli Half Wave Dipole In Free Space vs 6 cm Monopole Over 


Ground Plane 

Near field plots for magnitude and phase for the two 
antennas are presented in Figures 40 through 45. The near 
electric field is displayed in a.section of a plane in 3- 
dimensional space. For comparison purposes, the same view 
plane in the x and z directions are used for all such plots 
in this study. The origin of the view space plane is very 
close to the wire antenna and box surface (just 0.054 away) 
and extends to a distance of 10, gin the x and z directions. 
For these two simple antennas, additional view planes were 
selected in order to observe how the E-field varies away from 
the antenna. For this case, the results for peak E-field for 
the dipole are presented in Figures 40(a) and 40(b) (0.5, and 
0.8A away) and for the monopole in Figures 43(a), 43(b), and 
43(c) (1A , 2A , and 8A away). The electric near fields are 
referenced to one volt per meter for all plots. 

It is well known that the A/2 dipole has the same 
behavior as a A/a monopole (note that we test a ys monopole 
in this study). 

Before comparing the dipole and monopole field plots, 


consider the following: 
2p 
It 1s Known that AN eta Assume resonance for the two 
2 


antennas, then: 
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C250) 


The gain of an antenna over an isotropic radiator can be 
expressed as: 


2 
S E? / 
Ca = = ele 
Pp/4mR2 Pp/4mR2 
where S = power density (watts/m2) 
and PR = power radiated (watts) 
} Ep*/7] 
Then for the d/2 dipole: Ch ae ee (2.4) 
Pp/4mR2 
Ey?/7 
and for the monopole: Cm C225) 
Pp/4mR4 
#e is well known that: Gy = 2Gp (226) 
(In practice Cimaee oe Band |p 5 5B?) 


Forming the ratio of Equations 2.4 and 2.5 and using Equation 


2.6, we get: 








Gp Ep 1 

— = 5 i ee LS (257) 
GM Ey = 

However Dp = Ip? Rp 


Py = Iy* Ry 


In order to compare the gains of the two antennas, we have to 


Pp = Py : Ip* Rp = Iy? Ry (27.8) 


have: 


ag 


Using Equation 2.3 in Equation 2-3, we gee 


Rp 
2 e 
Lp Rp = Inu Thule Thus its 
2 








Thus, the electric fields will be equal only if the currents 
are equal for the two antennas. 
In Table 14, characteristics of the two antennas from 
the NEC output files are presented. From Figures 40 and 43, 
we conclude that the shape of near electric field is the same 
for the dipole and monopole antenna. The ratio of the input 
maximum current for the two antennas is: 
0.9764 X 1072 
20 log3g ——————_ = -7..28 GB 
0.2259 x 1071 
From Figures 40 and 43, we see that the O dB contours for the 
dipole are at a distance 2A away on the x-axis and for the 
monopole at a distance of 3) on the same axis. It appears 
that the field attenuates more rapidly with radial distance 
away from the dipole as compared to the monopole. For both 
antennas, the maximum value appears at a 40° angle from the 
x-axls. As the observation point moves away from the antenna 
along the x-axis, there is a difference of 4.6 dB between the 
electric field values on the contour plots of the two 
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antennas. Also, note that as the observation point moves 
away from the box surface along the z-axis, the difference 
appears to be 5.6 dB. Compared to the theoretical difference 
of 7.28 dB that we expect, there is a 2 dB offset which is a 
result of comparing a M5 monopole instead of a M4 one to 
the Ay2 dipole. 

Bhempnase plots of the E, and EZ components of the 
electric fields for the A /2 Gipole and 0.06 m monopole in 
Figures 41, 42, 44, and 45 show that both antennas have a 


smooth spherical wavefront pattern. 
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TABLE 14 
RESULTS FROM NEC OUTPUTS 














6 CM MONOPOLE 
OVER PERFECT 
GROUND 
(5-SEGMENTS) 


6CM MONOPOLE 
MOUNTED ON 
SURFACE PATCH 
BOX CENTER 
(5-SEGMENTS) 


LENGTH 
(meters) O. 75m (7/72) O.06m (775) 0.06m ( 4/5) 
WAVELENGTH 
Oe om Oc sm 
Vole LaVvoke 


DIPOLE 4/2 
IN FREE 
SPACE 
(5-SEGMENTS) 












(meters) O23m 












EXCITATION 
(voltage 
Source) 







FEED 

POSITION Sra lst lst 
# OF 

SEGMENT) 







MAXIMUM 
CURRENT 0.98 X 1074 0:22.% 10s- 0.26 X 1¢me 
(magnitude 

in amps) 


IMPEDANCE 
(ohms) 89.61 + 349.57] 26 - 435.82 | 10.99 - 437.53 
Ou4eee Oe 0.66 X 1072 0.36 X 1074 
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Figure 40. Total E-Field Contours 
Dipole A/2 on Z Axis in Free Space 
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Figure 40(a). Total E-Field Contours 
Dipole A/2 in Free Space (0.5A away from the wire) 


84 





Peony eee (be ner TO LVM) 
peewee GlnenaeeN Z-Axis 
eee Ota cme | CZ 


0.1856 


0.1392 


Z-AXIS VARTATION (METERS) 
0.0928 


0.04164 





0.0000 


0.0150 0.0614 0.1078 0.1842 0.2006 0.2470 
X-AXIS VARIATION (METERS) 





Figure 40(b). Total E-Field Contours 
Dipole A/2 in Free Space (0.8 A away from the wire) 
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Figure 41. X-Component, E-Field Phase Contours 
Dipole A/2 in Free Space 
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Figure 42. Z-Component, E-Field Phase Contours 
Dipole A/2 in Free Space 
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Figure 43. Total E-Field Contours 
Monopole 6cm Over Perfect Ground 
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Figure 43(a). Total E-Field Contours 
Monopole 6cm Over Perfect Ground (14 away from the wire) 
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Figure 43(b). Total E-Field Contours 
Monopole 6cm Over Perfect Ground (2A away from the wire) 
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Figure 43(c). Total E-Field Contours 
Monopole 6cm Over Perfect Ground (8A away from the wire) 
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Figure 44. X-Component, E-Field Phase Contours 
Monopole 6cm Over Perfect Ground 
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Figure 45. Z-Component, E-Field Phase Contours 
Monopole 6cm Over Perfect Goound 


93 


Ze Monopole at Center of Patch Box vs Monopole Over 
Ground Plane 


Input data files for NEC solutions of the near 
electric field of the monopole on the box center are 
contained in Appendix B. In Figures 46 to 50, the magnitude 
of the total electric field and the magnitude and phase of 
the E, and E, components are shown. In Figure 46, the E- 
field zero dB contour is observed to be closer to the antenna 
(it is at 2A vs. 3A) than for the monopole over a ground 
plane. Using the ratio of the maximum input currents alone, 
we calculate that the theoretical difference in the electric 
fields of the two antennas should be -1.08 dB. The contour 
plots of Figures 46 and 43 show that there is a difference of 
-3 dB along the box surface direction anda difference of 
+8 GAB up the z-axis, perpendicular to the surface, compared 
to the isolated monopole reference case. These differences 
in near field densities are attributed to box radiation and 
diffraction. In this near field region we can observe beams 
and nulls forming. A maximum occurs at 60° up from the box 
surface vs. 40° for the monopole over the ground plane. The 
main lobe (Figure 46) starts to develop at a distance 2) 
(0.6 m) from the antenna, along the x-axis. Beyond this 
point, it generally has the same pattern independent of 
distance. Within the 2A distance, the pattern shape has not 


yet fully developed. 
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In Figures 47 and 49, Ey, and E, contour plots 
demonstrate the behavior of the magnitude (the Ey component 
is negligible and is not presented). 

Recalling the Epeak field variations and noting that 
Epeak consists of E, and E, components, we can scan Figures 
47 and 49 and see how these components contribute to the 
total field of Figure 46. Different variations with distance 
along the x- and z-axes and with take-off (elevation) angle 
are noted. The z-component is predominant, as it should be, 
for a z-directed monopole. Null effects are most apparent 
for E, at moderate elevation angles and, as expected, E, 
nulls at the zenith. 

Figures 48 and 50 show the phase of the x and z 
components. The electric field phase has the same spherical 
wavefront pattern as the monopole over the ground plane 
(Figures 44 and 45), except for the region of the null, 
caused by the box. The disturbance in the electric field 
phase appears at the same 279 location in both phase plots 
(Figures 48 and 50). It is quite severe in the z component 
phase (Figure 50) and is barely observed in the x component 
phase (Figure 48). 

In order to get better insight into the electric near 
field variations and where the maxima and minima occur, 
3-dimensional plots are presented in Figures 51 and 52. These 
plots were produced from Fortran program III, described in 


Appendix B and use the same view plane as previous plots. 


95 


Each plot displays a surface whose elevation points 
represents field strength in the upper portion, and a normal 
2-D contour plot “projection” Gin the lower poneroen. Two 
different view points were selected, one "looking away from" 
the monopole and one towards it. In each figure the vertex 
that corresponds to the pointed "spike-like" area of the 
surface is the origin where the antenna is mounted. The 
delay of the field as the observation point moves away from 
the monopole source is expected. The null "trough" is very 


striking in this type of presentation. 
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Figure 46. Total E-Field Contours 
Monopole at Patch Box Center 
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Figure 47. X-Component, E-Field Contours 
Monopole at Patch Box Center 
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Figure 48. X-Component, E-Field Phase Contours 
Monopole at Patch Box Center 
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Figure 49. Z-Component, E-Field Contours 
Monopole at Patch Box Center 
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Figure 50. 
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Figure 51. Total E-Field 3-D Plot, View Toward Monopole 
(Monopole at Patch Box Center) 
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oe Monopole at Edge of Patch Box vs. Monopole at Center 


Figures 53 to 59, for the edge-mounted monopole 
follow the same format and trend as the center-mounted case. 
The following comments will amplify differences from the 
center case and interesting features of the field 
distributions (NEC input data sets are in Appendix B). 

The edge-mounted geometry provides a larger planar 
surface (box top) in front of the monopole view plane and 
results in the following differences with respect to the 
CencCer-Mount. Ggeenerry- 

- The z-axis peak field is somewhat greater at given 
distances from the box surface. This is traced to a 
larger E, component. 

- The elevation plane null is not as deep. 

- The phase contours for the z component shown evidence of 
two close-in nulls, but very little "phase wrinkles" 
beyond the near zone. 

- The contour plots depict a more uniform field overall, 


with a less severe null. 
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Figure 54. 
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Figure 57. 2Z-Component, E-Field Phase Contours 
Monopole at Patch Box Edge 
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Figure 58. Total E-Field 3-D Plot, View Toward Monopole 
(Monopole at Patch Box Edge) 


£10 


| 





yvev ® 
youl! 





(SSHOLWd LX 


| 


\)X 


s9q (Say: 
wa 


(Y4INIO NOMI WO £9°¢) 


Jas IOP Se GEM 


“IOdONOW 


Ne) (eS Ucaicnelel eM ened 


a ee ee re a 






Qala DYLOsT3 


Viewed From Monopole 


(Monopole at Patch Box Edge) 


Total E-Field 3-D Plot, 


ieeure O°. 


iiplok 


4. Monopole at Corner of Patch Box vs Monopole at Center 


The corner-mounted monopole near field plots are in 
Figures 60 to 66. The NEC input data is in Appendix B. The 
fields for this case fall in between the center and edge- 
mounted field configurations. That is: 
- The elevation plane null is between the other two 
Configuraeions” nuilic. 
- Phase wrinkles show one null, as does the center-mounted 
case, but it is not as severe. 
- Contour plots indicate a "flatter" field pattern with a 


fairly uniform distribution away from the monopole. 
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X-Component, E-Field Contours 
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Monopole at Patch Box Corner 


Figure 61. 


Ullb ee 


5 be" | 
(SQL IW NOLLUTMHA SEXU- Z 





cbuy U 


(METERS) 


X-AXIS VARIATION 


OS!1°0 


E-Field Phase Contours 


Monopole at Patch Box Corner 


X-Component, 


Figure 62. 
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Figure 64. Z-Component, E-Field Phase Contours 
Monopole at Patch Box Corner 
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Figure 66. 
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G. RADIATION PATTERNS 
Radiation patterns were calculated and are presented in 
Figures 67 and 68 (Monopole at center), 69 - 71, (Monopole 
at edge), and 72 - 74 (Monopole at corner). 
1. Monopole at Center 
In Figure 67, the vertical pattern shows that the 
maximum gain is very close to 5.15 dBi, the theoretical value 
for a monopole over an infinite perfectly conducting ground 
plane. The metallic surface of the box which re-radiates 
affects the pattern causing four lobes to appear instead of 
the usual two. In Figure 68, the horizontal pattern shows 
omnidirectional radiation from the box-monopole, which is 
fairly small electrically, and not expected to contribute 
much directionality in azimuth. 
2. Monopole at Edge 
The vertical pattern (x-axis cut) in Figure 69 shows 
an unsymmetrical pattern with the maximum lobe of 6.6 dBi 
skewed away from the monopole. The y-axis cut in Figure 70 
shows an expected symmetric vertical plane pattern. The 
horizontal pattern of Figure 71 shows a skewed pattern, with 
the maximum field opposite to the monopole location. 
3. Monopole at Corner 
The results of vertical and horizontal patterns for 
a corner-mounted monopole (Figures 72 - 74) are similar to 
the ones for the edge-mounted case with unsymmetrical 


patterns occurring. 
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Figure 67. Vertical Pattern, Monopole at Patch Box Center 
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Figure 69. Vertical Pattern (X-Axis Cut) 
Monopole at Patch Box Edge 
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Figure 70. Vertical Pattern (Y-Axis Cut) 
Monopole at Patch Box Edge 
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Figure 72. Vertical Pattern (X-Axis Cut) 
Monopole at Patch Box Corner 
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Figure 73. Vertical Pattern (45° Cut) 
Monopole at Patch Box Corner 
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meee WIRE GRID MODELING FOR NEAR ELECTRIC FIELD PREDICTIONS 


A. GUIDELINES FOR WIRE GRID MODELING USING NEC 

Solid surfaces can be modeled in NEC with a grid of 
wires, with the restriction that the grid cells are to be 
small in terms of a wavelength. Wire grid modeling 


guidelines are given in Reference 1. 


B. WIRE-GRID BOX MODEL FOR NEAR ELECTRIC FIELD INVESTIGATION 

For the wire grid modeling technique, the NPG1000 version 
of NEC allows up to 1000 wire segments. Typical run-times 
for wire grid models of the box-monopole are six times those 
of the surface patch models. The box of Figure 3 is modeled 
as a five-sided wire-grid box of 0.1m ( 3 at 1 GHz) per 
Side and cells of 0.0125 by 0.0125 meters. The 0.06 m 
monopole antenna was divided into 5 segments and placed on 
top of the wire grid box at the center, edge (3.75 cm from 
center), and corner (5.3 cm on the diagonal). The geometry 
of the wire grid box is shown in Figures 75 - 77. Appendix C 
contains NEC input data files for these geometries. The 
antenna was fed at the base segment for all cases. The wire 
grid geometry of Reference 8, was used for calculations of 
near electric field in the present study. This geometry 
produced good results for admittance and average gain 
compared with experimental data of Reference 5. The Fortran 


programs of Chapter II are also used for magnitude and phase 
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contour plots, and 3-D plots of near electric field former. 
wire grid box case with the same field point locations used 


in the surface patch model. 
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Figure 76. 
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MONOPOLE 6 CM AT CORNER (5.3 CM) OF WIRE GRID BOX 
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Monopole at Corner of Wire Grid Box 


Figure 77. 
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1. Monopole at Center of Wire Grid Box vs. Monopole at 


Center of Surface Patch Box 

The excellent agreement of near fields (Figures 78 to 
84) for the wire grid box with center-mounted monopole to 
those of the surface patch/center monopole model (Figures 46 
to 52) attests to the equivalence of the two numerical 
models. Differences are less than 1 dB, a value which is 


difficult to measure. 
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Figure 79. 
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2. Monopole at Edge of Wire Grid Box vs. Monopole at 
Edge of Surface Patch Box 


The edge-mounted monopole for the wire grid model 
compares well to the surface patch version. The differences 
are expected to be attributed to a more accurate rendition of 
edge effects for the grid model versus the patch model. 
Observable differences between Figures 85 to 91 and 53 to 59 


relate to deeper nulls for the wire grid model. 


142 


cs ~ . x ‘ ‘ ‘ ‘ 
es ‘ N ‘ aU . N ‘ ‘ 1 
a 
“S mn ‘ ‘ Sia Cs ‘ ‘ ; 
~ . . 
~ NN Rs oes Ne Sae¢ ‘ A 1 De 
Se = ‘ < Cn nr op Se ‘ ; a 
nin ‘ S x bee. ° 4 ‘ ; ae 
© x . ‘ A 
oR bs BS * ‘ S * H _ ‘ an 
x 
< Od, : . ‘\ ue ae i ‘ s Y 
~ 
7% eee ’ 
aS Be Bey A : N any 
fea) oe O-» ia ee ‘ ‘ © 1 U 
~—= ve a ane Sc aS “ . NY x wr ; “ 
= f . aS N M ‘ 
Ch . ce Ne ae cae . ‘ 
NOs ~ . Sy OS ’ ‘ ‘ ’ 
1.) _ . Or Ss Seo 8 ‘ ct 
aes . ~ 
‘ts & N ' (an) 
~ . Ny cS j te ‘ - 
f. ~ S SS Ss Si oe \ Pe wn 
GUEZs Pie. igs ‘ ess come . 5 ‘ 
(J MS ” S - \ ‘ : (an) 
. Bin S = \ ‘ © 
ae) PS Ny Sees a \ \ tm 4 
. re SeR yo N ‘ at { 
wy TS s oS ae aS a ? 
td . ~ ‘ By A 
a= @) are Sg SN Ne CPA iG i : 77 
ac é a se nee x Ny Sa ‘ ; oa 
TRE i ors x Ne 4 LY 
us ese ae << Ni etek Rw, ‘ G 
Se J - cairo ‘N x a) rane * ‘ if 
ae sae Se . 
as g ie a ae 
- * 
e151) Vor ay, SSN ees , t 
a ae ee 4 GY ' ‘ 
Cow < ~ ~ ~ LS ‘ ‘ 
GIN aa eS AMES vw ' 
a . ON 
. ae Py x N Ly a ‘ 
Ze A x“. Ss ay S ‘ J 
t-~fy) Be ~ mh % i’ \ A ‘ 
Ja -- ‘ . Ss, sos vot ‘ 
s = Ob Ol a x Se a ‘ ; 
4 4 oie Sate SR Deed 76 en ee Y ' ’ 
[1 J2.™m ae Bat PI ree ie SO 0 5 t 
Cit ; BA So eee ie ee 
. ‘ 
1409 of Se Se ay eo 
_ , oe = ‘ Gir, 8 4 
BF 
1t AG eee ee ee ; 
1c) ¥ ne Pees ISS Oo 4 
: wai - ae Sse ry 
It . ys a & Sect 
- ¢ ~ ie 1 
af : a ~ Se ee 
( a v,? ae OY 
i. . 
IN a) “1 poss x ' 1 
° ~ ~ ' 
9S: © — je i eer 
\ 7 ‘ Sus u 
7 ¢ --7 
a ras “ [0 +o 
te ‘ 
* ¢ 
é 
moe go 4 — 
ca 
¢ oa 
¢ 
5O" ne Pt 
ne : ’ 
- 
ra i , a“ 


Slo°e SLb°? Gon | G40" | S09" 0 SiO 
(SMW NOTLURIWA SIXU: Z 


2.679] 


2.299) 


TERS) 


~ 
r 


Le ig 


143 


Teles) 


X-AXIS VARIATION (M 
Total E-Field Contours 


Monopole at Wire Grid Box Edge 


Figure 85. 








O98b ‘¢c 


Bilb'e 


OAL I 
(SUL IW) 


ate! 


NOULUTSUA STXU- Z, 


a) S, 


yy 
Ld 
N 
t 
4 
‘ 
N 
t 
' 
4 


cbus 


0 





X-AXIS VARIATION (METERS) 


E-Field Contours 


X-Component, 
Monopole at Wire Grid Box Edge 


Figure 86. 


144 





Sete Via Ni fe t C= a" —. Kl ee to m4 — Pam —_— (ot Se ta 
SLASE OF YX COMPONENT OF S-FIE: 
' | Re CNT fas Sal — - a ry t ee ed ed Fe 
font rane ) 
Boe IK 
Ing 


2.9860 


{ 
{ 1 t 
SNS eee Beye Pte. (Or Wilkes Eo) 
MIM te leroy. Re Boe Aye 1am A 
= 


2.4418 


_ 
- 
w 
nn 

' ont 
ces be 
tae 
age Go 

\ ~- 
Gre art 
— 

1 — 
_— 


Z-AXIS VARTAT TON 
1.2634 


0.6892 


0.1150 





1.1275 1.7017 2.2759 
X-AXIS VARIATION (METERS) 





Figure 87. X-Component, E-~Field Phase Contours 
Monopole at Wire Grid Box Edge 


145 


i 

22m 
ae) 
--~ 09 


ont 


{a} ht 


Me. oSc 
ee Be od 
S kentbention 

ON cme tet 





OYUL 2 


BEitb dc 


ay 
Uocldeind 


by UCT 
NOPTUTWA SEXO: Z 


cbuy 0 


—— 


OSII'O 





2.850) 


ie7Cl7 


X-AXIS VARIRTION (METERS? 


1.1275 


E-Field Contours 


Monopole at Wire Grid Box Edge 


Z-Component, 


Figure 88. 


146 





Ole ogub*Z wlib'2 GL. | be9C" | Zbuy 0S11'0 
al (SMOLWL NOTLULYUA SIXU Z 





7037 2.2759 2.8501 


1.1275 
ah ales iea (Ce 


TERS) 


-— 
lod 
= 





E-Field Phase Contours 


Monopole at Wire Grid Box Edge 


Z-Component, 


Figure 89. 


147 


(Oy 


=D ELECTRIC FELD (08 RE 
MONOPO 5 
SE ON RecN WISE GID BOX 





aad WYLIT 


Figunes30 Total E-Field 3 
; iel D Plot, View 
) towar 
(Monopole at Wire Grid Box Edge) d Monopole 


148 


(MSINIO WONT WO SZC) 


XO8 CeO SYM 3O JO0d4 IV WO9 4 IOGdONOW 
NO s@)RGRE ERC ls 


Q13l4 HYLI313 





Viewed From Monopole 


(Monopole at Wire Grid Box Edge) 


Total E-Field 3-D Plot, 


Figure 91. 


149 


3. Monopole at Corner of Wire Grid Box vs. Monopole at 
Corner of Surface Patch Box 


The previous two correlations between surface patch 
box and wire grid models continue for the corner-located 
monopole. Values of 1 to 1.5 dB in field strength 
differences, with deeper nulls, summarize the results of 


comparing Figures 92 to 98 with Figures 60 to 66. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The goal of this investigation was to accurately predict 
near fields using the Numerical Electromagnetics Code (NEC) 
on ship-like structures. The near fields of Navy HF antennas 
are primary contributors to HF radiation hazards (RADHAZ). 
Providing reliable interference-free communications is a 
formidable task because of the high density of 
electromagnetic sources on the topside of a modern warship. 

Since no validation benchmark results for near fields 
were available for ship antenna installations, a modeling 
exercise was undertaken using NEC in order to select the 
optimum models which would be used to calculate near fields. 
Box-like structures were analyzed using surface patch and 
wire grid modeling techniques. The simulation models consist 
of a A/5 monopole mounted on the top of a N73 box at three 
different locations: center, edge and corner. Optimum models 
were selected by varying the patch density on the top surface 
of the box and observing the convergence of the solution for 
input admittance and average power gain, in order to ensure 
the validity of the models and improve the confidence in near 
field predictions. Optimum NEC models for the three 
different mounting geometries are: 

- CENTER: 7x7 = 49 patches on top Goole tepals area) 
- EDGE AND CORNER: 11x11 = 121 patches on top (.0009)2 
patch area) 


1S 


NEC results for input admittance correlate very well with 
measurements. Average power gain varied between 1.80 to 2.02 
(within the 10% engineering tolerance of the theoretical 
value). Single Precision (NPGNEC) calculations gave the same 
results as Double Precision (DNPGNEC). 

Even though edges are not modeled in the surface patch 
technique [Ref. 1], this study proves that, for positions 
very close to an edge, good results can be obtained by 
careful subdividing (no special subdivision of smaller 
patches in the vicinity of the edge or corner was required). 

Fortran algorithms were developed to produce near 
electric field (magnitude and phase) contours and 3-D plots 
using DISSPLA [Ref. 9]. The near field for the monopole on 
the box has the same characteristics in magnitude and phase 
as the monopole over a ground plane except in a region where 
nulls occur from box radiation and diffraction effects. The 
edge/corner-mounted geometries produce slightly different 
near electric field maps compared to center-mounted geometry. 

Surface patch and wire grid models for NEC gave similar 
results for near fields. 

Previously, generalized guidelines for near field 
modeling had not been developed for NEC and the use of wire 
grid and surface patch modeling for near field parameters was 
approached with caution. Guidelines developed in this study, 


as well as the results of the near field behavior of the 
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monopole antenna on the conducting box, can be used for 
future investigations on "sShip-like" structures. 

Recommendations for continuation of the present study 
are: 

- Extend the model to larger and more complex structures 
which more closely approximate the ship's topside super- 
StEbucCeune. 

- Develop additional Fortran algorithms for use in near 
field investigations. 

- Exercise a new version of the PATCH Code, "JUNCTION" 
[Ref. 11] in order to check the validity of the code 
and the computer/time cost in comparison to NEC. 

The present study is an important step in the direction 
of modeling the effects of shipboard environments on the near 
field of antenna structures. It demonstrates how NEC surface 
patch and wire grid modeling techniques can be applied for 
accurate near field predictions of antennas in complex 


environments. 
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APPENDIX A 
NEC INPUT DATA FILES 


PART A: NEC INPUT DATA FILES FOR 
SURFACE PATCH MODELING 


Pile. Seacnuo0 DATA Bl 


Cii SP30X 5X5 PATCHES ON TOP SURFACE 
ColorPecwwienc. | 00) GHz 
Sito oy Ue y= 0. 55 9,025,—0.551 


Sco, 90,0, -0. 5, 1 

Sispoe ls osu, 0, 045,050 
SHOU hi) eohe b(t ak 

GX 0,110 

Sloss +o oy 0551, 02.5, —0 Soak 
Seer peo ooo 1 

CS 0505051 

Geel 

Ee 0703050, L000 

Ci: 

WG 

xg 0 

Bel 


1. Data File Surface-Patch Box Structure 
(5 x 5 patches = 25 patches on top 


Surtiace). 


PEee Sr Sext DATA M1 


Cia SPBOX 7X7 PATCHES ON TOP SURFACE 
Ce PREQUENCY "1.010 GHZ 

215, 670.597 -055,0,0.,5,=0,5, 1 
360;0;,0,—-025, 1 

Silo; 550.0; 70.5700. 5, 0,0 
2€0,0,0.5,0.1 

GX dU, 110 

SU Gin =O), 5) 05. 0 Si, Oe Set 
scl07d,0.5,0.5,2 

Gs 0,0,0.1 

Ge 1 

FR 0,0,0,0,1000 

GN l 


2. Data File Surface=—Patch Box Structure 
(7 x 7 patches = 49 patches on top 


surface). 


alist 


PILE: oF bl DATA M1 


CM SPBOX 9X9 PATCHES ON TOP SURFACE 
CE FREQUENCY 1.00 GHZ 
SM355,0.5,—-0.5,0,0:5,-0 e508 
S€03070--0 0551 

SMS. 550.5. -05e0700 50h 0 
$c0,050.5,0,1 

xX Q, 110 

SM9, SP} 5 10) 6 5,70. Sele es 5,-0.4 5, a 
Se60,0,0¢5, 0.504 
GS 0,0,0.1 
Ged 

FR 0,0,0,0,1000 
GH 1 
HG 
XQ 0 

EN 


3. Data File Surface-Patch Box Structure 
(9 x 9 = 81 patches on top surface). 


FILE: SPBOX11 DATA M1 


CM 6 Ci MONOPOLE AT reduces OF SURFACE PATCH FREQ. 1.00 GHZ 


CM TOP SURFACE 7X7 PATCHE 
CE CALCULATION OF AVERAGE > GAIN AND INPUT IMPEDANCE 


GF 
GH1,5,0,0,1,0,0,1.6,0.016 


GS 0,0,0. 1 

GEL 

Xx 0,1,1695159 

RP 0,31,4,1002,0,0,35,15 
XQ 0 

EN 


4. Data File 6cm Monopole at Center of 
SPBOX. Input Admittance and Average 
Power Gain Caliculations 
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FIEE: SPEDGE] DATA me 


Sor eOcc. stl le rAnTecHesmOnN TOF SURFACE (6CM Al EDGE 3.65CHM FROM CENTER) 
Ge FREQUENCY 1.00 GHZ 

SHSn Sec Ss Soe Sra, Spl, Bis at 

BeUs0,0,-0.55 1 

SHS 5 Sp Ue So —Ws Se Bla) oy alles) 

eeu,0,0.5,051 


ox 0,110 

Sil, 2, =0.5,-0.5,1,0.5,-0.5,1 
5cO,0,0.5,0.5, 1 

Gs, 0,0,0.1 

SE 1 

FR 0,0,0,0,1000 

GH 1 

WG 

XQ 0 

EN 


5. Data File Surface-Patch Box Structure 
Cis ti] 221i patches on top ustiriace ). 


Bebe: SPEDGE]1 DATA M1 


cr 6 Cm MONOPOLE AT EDGE 3.63CM FROM CENTER FREQUENCY 1.00 GHZ 
Cr POP SURFACE 11X11] PATCHES 

CE CALCULATION OF AVERAGE GAIN AND INPUT IMPEDANCE 

GF 

GW1,5,-0.363636,0,1,-0.363636,0,1.6,0.016 


GS 0} 510) 5 (0), 3 

Geol 

x O;l>1,0,1,0 

RP 0, 351,13,1002,0,0,3,15 
XQ 0 

EN 


6. Data File 6cm Monopole at Edge (3.63cm from 
Center) of SPBOX. Input Admittance and Average 
- Power Gain Calculations. 
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FILE: SPCNi Rieti M1 


CM SPCNTR 11X11 PATCHES ON TOP SURFACE (C6CM MONOPOLE AT CORNER 5.14 CM 


CM ON THE DIAGONAL FROM CENTER 
CE PREQUENCY S00 Giz 
SM3.5,0. 55-02 5500.50 oe 
SC0n0e 0, — Ceo 

SMS 55,0 .5s Uso 70> Ure. Un 
SCO; 070 25),.0\77 

GX 0,110 

SMl1 3115-075. -025, 15025, 00 50 
$C05070.52 02572 

Gs 0510, 01 

GE) 

FR 0-0,0,0, 1/000 

GN ol 

iG 

XQ 0 

EN 


7. Data File Surface-Patch Box Structure for 
Monopole at Corner of SPBOX. 


FILE: SPCNT RIL DATA M1 


CM 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
CHM TCP SURFACE TIX) 1 PaAneaHes 


CM FREQUENCY 92-00 | 3GHZ 
CE CALCULATION OF AVERAGE GAIN AND INPUT IMPEDANCE 


GF 
GW1,5,0.363636,0.363636,1,0.363636,0.363636,1.6,0.016 


GS 0,0,0.1 

GE 1 

EX 0,1,1,0,1,0 

RP 0,31,13,1002,0,45,3,15 
XQ 0 

EN 


8. Data File 6cm Monopole at Corner (5.14cm on the 


Diagonal) of SPBOX. Input Admittance and 
Average Power Gain Calculations. 
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PART B: NEC INPUT DATA FILES FOR 
NEAR ELECTRIC FIELD IN 2-D 


Prove: SNGr DATA Bl 


Ci SPS8Ox7 x7 “PATCHES ON TGP SURFACE 
GE FREQUENCY 1.00 GHZ 

Sits oe SU On. Os MO ok 
SG0r. 0.0 Ole See 

Silay on. 5s =O. S09 0. 5,050 
$CO70,0.5,0,2 

Gx 0,129 

Sitio 75S eo edo, lO. 5,-0.5, 1 
SC 0p Oy 0. oc sk 

Gs 0505001 

Ger 

FR 0,0,0,0,1000 

GH 1 

WG 


NX 
C46 CM MOHNOPOLE AT CENTER OF PATCH BOX 


GE HEAR E=- LEU GN X-AXIS A? 0, 0166,0,0. 215 


GF 

(hh Ss pA le Wn Ola 2h on lo (lake 

GS 0 

GE 1 

EX 0 

Elec 
0 
0 


HMrHO fF 


1,0 
NE »1,0.0166,0,0.115,0.01,0,0 


XQ 
EN 


1. Data File Near Electric Field Along 
xX-Axis for Monopole at Center of SPBOX. 


FICE: SHERI DATA BL 


CM SPBOX 7X7 PATCHES ON TOP SURFACE 
CE FREQUENCY 1.00 GHZ 

Silos o> Dao 0D, Up 0. 5,-0.5) 

SE pW al, Slo Sip dt 

Suse SaWlcbn essen Wp Wes, 0, 

Seo, 0,0. 57 Uist 

GX 0,110 

SM7 a7 p70 25, 05D, Lp oe 0 wo, 

SC, Oy 0 saa orl 

GS 050,0.1 

GE 1 

FR 0,0,0,0,1000 

GN 1 

WG 

NX 

CM 6 CM MONOPOLE AT CENTER OF PATCH BOX 


ae NEAR E-FIELD ON Z-AXIS AT 0,0,0.175 
G 
GW1,5,0,0,1,0,0,1.6,0.016 


rue oOo 
e 


71,0 
0,0 


MNF © 


,0,0.175,0,0,0.01 


Z. Data File Near Blectric Field Along 
Z-Axis for Monopole at Center of SPBOX. 
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FILE: NEDGEX DATA Bl 


CM SPEDGE] 11X11 PATCHES ON TOP SURFACE (6CM AT EDGE 3.63CM FROM CENTER) 
CE FREQUENCY 1.00 GHZ 

SM3,35,0.5,-0.5,0-0.55—-0.551 

$60; 0, 05-05. 1 

SM3,3,0.5,-0.5,0;0.5,0,0 

SsC0;070. 5,054 

Gx 0,110 

SM11,11,-0. 5.70.57), 0255-0. 

SCO, 0, 0crGeo. ) 

GS De 


FR 0,0,0,0,1000 
1 


CM 6 CM MONOPOLE AT EDGE 3.63 CM FROM CENTER (11X11 PATCHES ON TOP) 
CE NEAR E-FIELD ON X-AXIS 


Gh 1, ee 363636,0,1,-0.363636,0,1.6,0.016 


,0, 
il 
2 
2 -0.0197636,0,0.115,0.01,0,0 


3. Data File Near Electric Field Along x-Axis for 
Monopole at Edge of SPBOX. 


FILE: -NEDGEN DATA Bl 


CM SPEDGE] 11X11 PATCHES ON TOP SURFACE C6CM AT EDGE 3.65CM FROM CEN 
CE FREQUENCY 2}. 00.GHZ 

SMS5 S70 oa - Oo 0, Ola 5S Oreo 

SCG, 0,0;=0.5,1 

SM575 0s op 05s 005 00 

SCO, 050.5, 0,1 

GX 0,110 

SM1LI, i), =0.5,-0.5,1,0.5; -0. 5-2 

sC0,0,055,0, 5,1 

Gs 0,0,0.1 


FR 0,0,0,0,1000 


QnG 
=70m 
Hor 


CM 6 CM MONOPOLE AT EDGE 3.63 CM FROM CENTER (11X11 PATCHES ON TOP) 
CE NEAR E-FIELD ON Y-AXIS 


GH 1, Se 5636356,0,1,-0.365656,0;51.6,0.016 


20. 
uh 
5 
2 


Ne ot 


-0.0363636,0.0166,0.115,0,0.01,0 


4. Data File Near Electric Field Along y-Axis for 
Monopole at Edge of SPBOX. 
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FILE: NEDGEZ DATA Bl 


CM SPEDGE1] 11X11 PATCHES ON TOP SURFACE (6CM AT EDGE 3.63CM FROM CENTER) 
GCE FREQUENCY 1.010 GHZ 
Bho so5 0.5, —-0.5570,0.5,70.5,1 
se0,0,0,-0.5,1] 
prov o, 0.5,-0.5,0,0.5,0,0 
5c0;,0,0.5,0,1 
Gx 0,110 
Sele tl, =0.5,-0.5,1,0-5,-0.5,1 
ee0,0,9.5,0.5,1 
GS 0,0,0.1 
GE 1 
FR 0,0,0,0,1000 

1 


CM 6 CM MONOPOLE AT EDGE 3.63 CM FROM CENTER (11X11 PATCHES ON TOP) 
CE NEAR E-FIELD ON Z-AXIS 


GW 1,5,-0.363636,0,1,-0.363636,0,1.6,0.016 


, ,0 


rO,O 2 
ple Lei 

pay Ire 
»1,1,200,-0.0363636,0,0.175,0,0,0.901 


v 
rm 
Se 


5. Data File Near Electric Field Along z-Axis for 
Monopole at Edge of SPBOX. 


EELE: NCORNX DATA Bl 


CM SPCNTR 11X11 PATCHES ON TOP SURFACE (6CM MONOPOLE AT CORNER 5.14 CM 
CM ON THE DIAGONAL FROM CENTER 
BES FREQUENGY 1.00 GHZ 
Sio,o,0.5,-0.5,0,0.5,-0.5,1 
5c0,0,0,-0.5,1 

SMopis Ole > sO. 5,0, 0.50, 0 

peor 0,0.5,0,1 

Gx 0,110 

Polls -0.5,-0,5,1,0.5,-0.5,1 
eeo,d,0.5,0.5,1 

Gs 0,0,0.1] 

GE 1 

FR 0,0,0,0,1000 


©) 
=< 
~ 


CM 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
CMOTOP SURFACE 11X11 PATCHES 

Ci FREQUENCY 1.00 GHZ 

ae NEAR E-FIELD ON X~AXIS 

G 

GW1, a “F Bee Pe 702563656. b Oe SO a cece aaa 


NN o- 
HHO ~ 


70.0529656,0.0565656,0.115,0.01,070 


6. Data File Near Electric Field Along x-Axis for 
Monopole at Corner of SPBOX. 


ow 


FIEE= NCORKT DATA Bl 


CM SPCHTR 111X111 PATCHES GN TOP SURFACE C6CM MONOPOLE Al CORNES SS ucees 
CM ON THE DIAGONAL FROM CENTER 

CE FREQUENCY 1.00 GHz 

Si1Si5 55 Ue SUS UU oo 0 ee 

SCO; 0,0, -0-5.. 

SMSe5> 0.55 -On oo U0. 5) UO 

SEO) 5 ll 5 18) 5 Sa (ls 2h 


Gx 0,110 
’ 05520555 150255 -02 551 
065.4 


CM 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
Cit TOR SURPACE ExT) (PAleGnes 

CM FREQUENCY 1.00 GHz 

CE NEAR E=FIELD On yoexts 


GF 
GH1L,5,0.363636,0.363636,1,0.363636,0.362636,1.6,0.016 


, , 


0: 
vt 
275 
1,2 


. 03656356, 0'.0529656,0.115,0,70-0270 


7. Data File Near Electric Field Along y-Axis for 
Monopole at Corner of SPBOX. 


FREE: NEORNZ DATA Bl 


CM SPCNTR 11X11 PATCHES ON TOP SURFACE (C6CM MONOPOLE AT CORNER 5. 1amer 
CM ON THE DIAGONAL FROM CENTER 
CE FREQUENCY 1°00 GHZ 

SMS, ors oe 0. Oy Oo, 0 
SG050- 0). —0r- 57 

SiS5 55. Ue oo Ul oer Oo Ue 
SC0p0,0.5,0;1 

Gx oO, 109 

Sl Sl is 0.5, -0.57 1-0 75,-0- 57 
S60, 050.5055) 

GS 0,090,001 

GE 1 

FR 0,0,0,0,1000 

GN l 

WG 


NX 
CM 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 


CM TOP SURFACE I1X1) PATCHES 
CM FREQUENCY 1.00 GHZ 
ne NEAR E-FIELD GN Z-AXIS 


G 
GHI1,5,0.363636,0.363636,1,0.363636,0.363636,1.6,0.016 


? 


a 


0 
00,0.0363636,0.0363636,0.175,0,0,0.01 


’ 


0, 0.1 
Deis 0 
2ooeo 
laloe 


8. Data File Near Electric Field Along z-Axis for 
Monopole at Corner of SPBOX. 
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APPENDIX B 
FORTRAN PROGRAMS AND NEC INPUT. DATA. FILES 


PART A: DESCRIPTION OF FORTRAN PROGRAMS FOR INTERFACING 
NEC OUTPUT WITH DISSPLA GRAPHICS PACKAGE 


PROGRAM I: CON1 

This Fortran program is constructed in order to generate 
contour plots of magnitude Epeax (Etotal) of the near 
electric field using NEC output in conjunction with the 
DISSPLA graphics package [Ref. 9]. Part. or He's typical 
output which shows the calculated Epeax, is given in the 
example 1 data file, and is produced using the PL card with 
@peions 2, 2, 5, 1 and an NE card [Ref. 1}. In example 1, 
the 2nd column is the x-variation, with points spaced as is 
defined by the NE card, while the third column is the 
magnitude of Eneak (Etotal) which corresponds to each point 
in space. 

The CON1 program reads the values from the 3rd column of 
NEC near field output, i.e., example 1 data file. It then 
arranges these values in an NX by NY array, for input to the 
subroutine CONTOR, which in turn creates the contour plots. 
CON1l also selects the maximum (Emax) and minimum (Emin) 
values of the electric field (3rd column of example 1 data 
le) . Using these values in combination with the NLEVEL 
value (the number of levels between the contours, which in 
the present study is twenty), it calculates the variable FINC 
which is the contouring interval, which is needed by CONTOR. 
The values of Emin, Emax, and FINC are displayed on the 


screen. The calling statement is: 
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CALL CONTOR (A, NX, NY, FINC) 


where: A - NX by NY array containing the values to 
be contoured 


NX - Number of columns in the array (in the present 
study 1 is 30) 


NY - Number of rows in the array (in the present 
study it 1s 30) 


The values in the array A have been arranged by CON1 as 
follows: 


Z A(1,NY) TGA AINBO Ns 5 aPPe Wap Gans) 
AG yo) 


5 IIS, AL) 





In the contour plots of study, the magnitude of Eveak? 
the electric field, is presented in a form of a section of a 
plane in 3-dimensional space composed of a (30 x 30) 900- 
point-window in space. The unit of magnitude used to express 
the electric near field in all the plots is dB with reference 
to one volt per meter: 


|E | 
adB(1V/m) = 20logi9 mae) 
1V/m 


When the magnitude of Epeak (Etotal) is plotted using NEC 
output in the form of example 1, the CON1 FORMAT is defined 
SS ee 0 3) ee However, when Epneak Of the inividual 


components (Ex or Ez) is plotted, the PL card with options, 


170 


oot, 1S used in Conjunction with an NE card [Ref. 1]. In 
Mars case, NEC output has the form of example 2. Columns in 
mie example 2 data file are: xX-variation, Ex, x-component 
Phase, Ey, y-component phase, Ez, z-component phase. In this 
case, when the magnitude of individual components, x- 
component (Ex) or zZ-component (Ez), is plotted, the FORMAT in 


Beomemust be: (15x, D1l0.3) and (55x, D10.3) correspondingly. 


CONTOR: This subroutine is constructed using the DISSPLA 
Graphics package [Ref. 9] and is modified for this thesis for 
near electric field contours. The subroutine has comments in 
each step explaining its function. The variables XORIG, 
XSTP, XMAX, YORIG, YSTP, and YMAX, have to be defined 
manually by the user for each case. POrmewcach plot, the 
program must be compiled by the "FORTVS" command before 
execution. In the present study, the parameters labeled 
y-axis are actually the z-axis parameters. 

The subroutine MYCON represents negative dB contours by 
dashed curves while the positive ones are solid curves. The 
humber in parentheses of the statements, COMMON WORK and CALL 
BCOMON [Ref. 9], must be 5,000 for magnitude plots while for 
phase plots it must be 20,000 because DISSPLA needs 


additional workspace for smoothing the phase values. 


PROGRAM II: CON3 
This Fortran program generates contour plots of the 


phase of individual components (x or z2 components) of the 


Ib fal 


electric field from NEC output and the DZISSPIA Grapes 
package [Ref. 9]. Part of a typical NEC output is presented 
in the example 2 data file. This output is produced from the 
the PL card with options 2,2,4,1 and an NE Card [Rete 
The columns presented in example 2 are: x-variation, Ex, x- 
component phase, Ey, y-component phase, EZ, z-component 
phase. 

The program reads the values from the 3rd and 7th 
columns (Phase of x and Zz components) and arranges them 
in the Array A consisting of Nx by Ny points using the 
FORMAT (25x, D10.3) and the FORMAT (65x, D10.3) for the phase 
of x and z~-components correspondingly. For the case of the 
phase plot in the present study, the array is selected as 100 
x 100 = 10,000 points. This is necessary because the near 
field region phase is very complicated, having abrupt 
variations and the DISSPLA program needs more points for the 
smoothing and accuracy. The number of levels between 
contours in the phase plots is selected to be: NLEVEL = 10 
versus 20 (which was used for the magnitude plot). The 
number in the statement: COMMON WORK and CALL BCOMON must be 
20,000, because in the case of phase plots, DISSPLA needs 


additional workspace for smoothing. 


PROGRAM III: CON2 
This Fortran program generates 3-D plots of magnitude 
Epveak (Etotal) of the near electric field. In this type of 


plot, a section of a plane in 3-dimensions consisting of 


LEZ 


(30 x 30 ) 900 points is presented. Electric field values 
are presented as a 3rd dimension (vertical axis). The 
contour plot appears at the bottom of the graph (exactly the 
same plot as the one produced by CON1) while on the top, the 
electric field variation as a 3-D surface is presented. The 
program is the same as CON1 except that the user must input 
the number of levels (NLEVEL) manually from the keyboard and 
also has the option to select coordinates of the viewpoint: 
wemev2Z, and VE. CON2 uses the subroutine DMH002 of DISSPLA 
(Ref. 9] which is listed at the end of the program. 
Parameters in the statements CALL GRAF3D, CALL GRAF and CALL 
RLVEC3 have to be input by the user and the program has to be 
compiled by "FORTVS" before execution. All other 
information, a user may need, appears in comments in the 
program. 

In order for DISSPLA to be automatically executed for 
each of the Programs I, II, or III, an EXEC file has been 
Penstructed correspondingly (CON1, 2, or 3 EXEC). The user 
has to input the filename of the NEC output data file (which 
has the form: FILENAME PLOTDATA) in the EXEC file and then 
typing CON1, CON2 OR CON3 in CMS, causes the execution of 
DISSPLA automatically. For all the programs the output is: 
DISSPLA METAFILE T, which can be copied to the A or B-disk 
and then, using the DISSPOP option, the user can select the 


type of output needed from a menu (Sherpa, IBM79, etc.). 


ES: 


FILE: DIPG 
N809 0.242E+01 
NedSa 0. 2526201 
N809 0.262E+01 
N309 0.272E+01 
N809 0.282E+01 
NS O92 eco On 
N809 0.166E-01 
N809 0.117E+00 
NoO9 022i 7 Er 00 
NS09 0.317E+00 
N3809 0.417E+900 
N809 0.517E+C0 
N809 0.617E+00 
N809 0.717£+00 
N809 0.817E+60 
N809 0.917E+00 
N809 0.102E+01 
NS OS ro, LPZEOL 
N809 0.122E+01 
N809 0.152E+01 
NS09 0.142£&+01 
N809 0.152E+0! 
N309 0.162E+01 
Néouo  OUs17Z2er01 
Heds “O-1éZ2E701 
N809 0.192E+01 
N809 0.292E+01 
N809 0.212E+01 
N809 0.222E+01 
N899 0.232E+01 
N809 6.242E+01 
N809 0.252E+01 
N809 9.262E+01 
N809 0.272E+01 
N809 0.282E+01 
N309 0.292E+01 
1. Example 1: 
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TFOMTDKDDODODOVDOODOO GOOD ODO DODD DODO DOOD O OOOO O 


PLOTDATA Bl 


SICSE—OT 
Jose —Un 
ao oe Or 
ooes0t 
-SHUE-01 
_2O3c= 0m 
sJLLE=02 
162GE=02 
el GeuF 
-L66E-01 
-oe E01 
ee ean 
pe yaN aS wal 
sso 9E=01 
-416E-01 
-461E=0T1 
-S04E=01 
-S45E—00 
-5384E-01 
alors) 110) 
6250 E=04 
-688E-01 
Pi SE=On 
67496E-01 
17 HE=0m 
of ISE-01 
-SL6OE=01 
-CoOE-Ut 
Oo SEOUL 
-oo0E=04 
soo SES UL 
-o95e>01 
TUDE —00 
-914E-01 
wIZec sur 
—t2OE [01 


NEC Output 


DIPH44 


0 

0 

0 

0 

C 

0 

0 

Cc 

0 

0 

0 

0 

U. 

0.144E+01 
0.147E+01 
Ue 50E+01 
Oel52E+01 
0 SLS55E+O] 
OF I5SSéE+Ol 
OF IGIEtOL 
0.164E+01 
U2167E+01 
0.170E+01 
Gels SEtOI 
O.176eE+01 
O77 et0l 
OSTS6LE+O1 
0.184E+01 
O21s7E+0) 
OF TIVE+O1 
O.19SE+01 
O.196E+01 
O2199E+0) 
GO. 202E101 
Decuoe +0] 
0.2082+01 
Oaclue+ol 
Neel SE+O! 
O.2l6E+ol 
0. 219E+Ol 
Og. 222e+01 
02 225E+01 
Orc2oetol 
UecslE+O1 
0.234E+01 
Oneoreroy 
O.259Er0] 
0.242E+01 
0.245E+01 
0.243E+01 
He2slE+01 
0.254E+01 
O2257E*01 
U22oUE+01 
OT 265E+01 
0.266E+01 
0.268E+01 
Oneal E Eo 
O,279E+01 
Ueer ETO! 
0.280E+01 
0.2835E+01 
).286E+01 
O,2s7Et+ol 


od 


NNNANNAHAAAAAN 
nu DRBWANRr OO 
PRONE ARAOUOWNO 


noo 9ES 
HO 2s 
267 Ses 


,6O0E= 


DNDNDDDAHDAADAAAAO 
OCwowwowowwowowowowowowvos 


OODOSOSDDGAGVOOOCVOGOCHSGVCABDDDOODCVOGOCCODOODODODODOODDOODODOOCOCOOOCOCOOCOCOOC COCO CO 
WCODWOrNWNHAMONANWOAin 


PLOTDATA Bl 


mmmmmm 
| 
ooo°o°eoo 
test 
ooo 


OOOO TDOCCDOCOCODDOCODDODOCDODCOCOCOOCOCOOFCOD OOO OMOOCOOOC OOO OOOO OOOO O 


ON 
154) 


ON 
+ ON 
N 


ON ON 
OOo 
in 

ih 


ON 
[ee] 


ON 
\o 
fo) 
di 
I 
OOGODDODDHOCOCOOCOCODOOODODOODOOCOD DODO OOOD ODO OOOOOOOO OOOO OOOOCOOO 
Len) i 4 


RNRAWNSN 
Oo OO CO 00 
Wino 


Bet tt et ee ee ee ee ee ee ee ee ee ep pe pp pp 


MMMM MMMMMMMMMMMMMMMMMmMMM 


-USO0E+0s 
Pl6 tao 
L7TGEFOS 
.174E+03 
-1l6]E0S 
.2149E+03 
plo SEs 0S 
mec E+ OS 
~108E+03 
~942E402 
PSVOEFOZ 
S65 5E+ 02 
~502E+02 
sooo erUe 
~2O5ErO2 
~505E+02 
PLO7ETOZ 
.266E+02 
MGAS Ech OZ 
2S93Et0Z 
-f S9ETOZ 
92SEt02 
SLTOEFOS 
~127ET03 
wa SELOS 
165EFO0S 
lV 9ELOS 
JelE+os 
s14SE+03 
»12SE+05 
.105E+03 
~865E402 
"67 LETO2 
247 7ET02 
peo LEtOZ 
.836E+01 
HlI6ETOZ 
~ SL8E+62 
-SLLETOZ 
sf e2OeE+02 
Wooo EtO2 
~114E+03 
. L135E+03 
.156E+03 
,I7SETOUS 
PLOLEFOS 
-ISIETOS 
~117E+03 
J 50E+O2 
-12eetoZ 
PS0SeFoc 
.280E+02 
~534E+01 
-176E+02 
~404E+02 
~635E+02 
~867E+02 
eELOETOS 
~1LSGETOS 
7157 ETOs 
ol7JETOS 
~LS5S5E+03 
~L31E+03 
PhO7 EOS 


2. Example 2: 


SOOSGOSOGDDSBDTVODDODOODODDODOCOCCOCOODDOOCODO OOD OOOO COODOCOOCCOC OOOO CO OOOO OOO 
eo) kel ie) fete) a o «© «© e jolene) © « «! «, («| °° 0« © «© © «+ «© @ ope ic sl ce e) <o “© of NeMMonMefiie) (sielteliise (s) :s) «) sis) wie mucine 


eS 


.000E+00 0.000£+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
O00E+00 0.000E+00 
OOOE+00 0.000E+00 
OOCE+00 0.000E+00 
QOOE+00 0.009E+00 
OOOE+00 0.000E+00 
OUO0E+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.00CGE+00 
OOOE+00 0.000E+00 
OUOE+00 0.000E+00 
OOOE+00 0.000E+00 
JOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOUE+00 0.000E+00 
OOOE+00 0.000E+00 
OO0E+00 0.000E+00 
OOOE+00 0.000E+00 
.OO0E+00 0.0C9E+00 
OOOE+90 0.000E+00 
.O0O0E+00 0.CO0E+00 
.OOO0E+00 0.000E+00 
.OOOE+00 0.000E+00 
OO0E+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOCE+00 0.000E+00 
OUCE+00 0.000E+00 
OOOE+09 0.000E+00 
OOOE+00 0.000E+00 
.OOOE+00 0.C00E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOCE+CO 0.000E+00 
OOOE+00 0.009E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OUOE+00 0.000E+00 
OOOE+00 0.000E+00 
.OOCE+00 0.000E+00 
.000E+00 0.0C0E+00 
.OO0E+00 0.000E+00 
QOOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
.Q00E+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
OOOE+00 0.000E+00 
QOOE+00 0.000E+00 

.000E+00 0.000E+00 

.000E+00 0.000E+00 
NEC Output 


OSOSOOSDSOSGCOCODODOOCODODODDODOODDOCOOCDOCODDOOOCCOOCOCO OC OOOOCOCOOCOCO OOOO OCOCOOOCOCCOCOCOe 
e e e e . e ° . ° . e . . . e . e e e e e . ry e e e e . e . ° . e . . ° e e e e e . e . . . ° . ° e e e e e e e e . * e e . ° . 


~460E+02 
PooGErTie 
~207E*+O02 
~769E+01 
-553E+01 
~190E+02 
sal A lea 02 
.466E+02 
~608E+02 
-TSLETO2 
~898Et02 
-LOSE+OS 
.120E+03 
~boSehOS 
PESUEtGS 
.166E+03 
~178E+03 
161E+03 
~145E+03 
~128E+03 
SLE + OS 
.943E+02 
oT TOUETO? 
.595E+02 
Pol 7 E02 
.238E+02 
-565E+0 
~27 E+ 02 
solSeEr0Z 
.500E+02 
~690E+02 
~881E+02 
-LO7Et0S 
-127E+03 
.1467E+03 
~167E+03 
-L7SE+05 
~TS53E+03 
~133E+03 
-L12E+03 
~ILZE+02 
. 7 02Z2E+02 
~491E+02 
~2/8E+02 
.635E+01 
~153E+02 
so LETOZ 
-590E+02 
~8l1E+02 
.1O3SE*03 
~126E+03 
.148E+03 
ole LEFOS 
.166E+03 
143E+403 
~L120E+03 
~968E+02 
.734E+02 
-498E+02 
e2oletue 
~-234E+01 
-2l6E+02 
-457E+02 
-699E+02 


FIEE: Cont FORTRAN Bl 
REAL¥8 RO 
C DEFINE THE ARRAY ACNX,NY) WHICH CONTAINS 
C THE VALUES TO BE CONTOURED.NX=NUMBER OF COLUMNS 
C AND NY=NUMBER OF ROWS OF THE ARRAY. 
C 
DIMENSION AC30,30) 
EMIN=1LE20 
EMAX=-1E20 


DO 50 I=1,30 
DO 40 J=1,30 
READ (€12,20)RO0 


C USE THIS FORMAT WHEN THE MAGNITUDE EPEAKCETOTAL) 
C IS 10 BE PLOTTED-NEC OCUTFUL IN als ease nao ane 
C BE PRODUCED BY PL 2,2,5,1 AND NE CARDS. 
Cc 
20 FORMAT (17X,D10.3) 
C USE THE FORMATS BELOW WHEN THE MAGNITUDE EPEAK 
C OF X AND 2 COMPONENTS CORRESPODINGLY 15 70 BE 
C PLOTTED.NEC OUTPUT EN THESEsGASEsanssen0eoe 
C PRODUCED BY PL 2,2,4%,1 AND NE CARDS. 
Cc 
C20 FORMAT (€15X,D10.3) 
C20 FORMAT (55x, D100 3) 


ACJ,1I)=20XLOG1OCRO) 
C IN THIS STEP THE PROGRAM PREDICTS THE MINIMUMS AND 
C AND MAXIMUM VALUES OF THE ELECTRIC FIELD. 


C 
EMIN=MINCEMIN,ACJ,1I)) 
EMAX=MAXCEMAX,ACJ,1)) 
G0 CONTINUE 
50 CONTINUE 
NX=30 
NY=30 
C DEFINE ue NUMBER OF LEVELS BETWEEN THe CONTOURS. 
NLEVEL=20 
C THE PROGRAM CALCULATES THE CONTOURING INTERVAL. 
FINC=CEMAX=EMIND7 CNLEVEE=ITD 
C EMIN,EMAX,FINC ARE DISSPLAYED ON THE SCREEN. 
WRITE €6,%) EMIN, EMAX, FINC 
C THE FOLLOWING STATEMENT IS THE CALLING SPATEMENT 
C OF THE SUBROUTINE CONTOR: 
CALL CONTOR CA,NX,NY,FINC) 
STUF 
END 


3. PROGRAM I: CON1 FORTRAN. Contour plots of 
Magnitude for Near Electric Field. 


FILE: Cont EXEC Bl 


FILEDEF 12 DISK FILENAME PLOTDATA A 
EXEC DISSPEA CONT 


4. CON1 EXEC: 
of DISSPLA. 


Fortran File for Automatic Execution 
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CON00010 
CON00020 
CON000350 
CON00040 
CON00050 
CON00060 
CON00070 
CON00030 
CON00090 
CON00100 
CONO0110 
CON00120 
CON00130 
CON00140 
CON00150 
CON00160 
CON00170 
CON00180 
CON00190 
CON00200 
CON00210 
CON00220 
ConNd0g230 
CONO00240 
CON00250 
CONO0260 
CON00270 
CON00280 
CON00290 
CON00300 
CON00310 
CON00320 
CON00330 
CON00340 
CON00350 
CON00360 
CON00370 
CON00380 
CON00390 
CON00400 
CON00410 
CON00420 
CON00430 
CON00440 


Bebe: CONTOR 


NQNNOAD ANANNANNAANANA AANNNANRO 


FORTRAN Bl 


SUBROUTINE CONTORCA,NX,NY,FINC) 


THIS SUBROUTINE CONTOURS AN NX BY NY ARRAY OF REGULARLY SPACED POINTS. 


A 
PSE 
NY: 


PING: 


AN NX BY NY ARRAY OF REGULARLY SPACED POINTS 
NUMBER OF POINTS Dh THE X=DIRECT ION 
NUMBER OF POINTS IH THE Y-=DIRECTION 


CONTOUR INTERVAL 


DIMENSION ACNX,NY) 
THE NUMBER IN THE PARENTHESIS OF THE FOLLOWING 
COMMAND DEFINES THE WORKSPACE IN THE DISSPLA 


GRAPHICS 


PACKAGE.THIS NUMBER HAS TO BE 5000 WHEN 


DHE MaAGHITUDE EPEAK IS 10 BE PLOTTED AND 20000 


WHEN THE 


PHASE OF THE INDIVIDUAL COMPONENTS IS 


MOSEE PLOTTED. IN THE CASe OF PHASE PLOT THRE ARRAY 
IS LARGER AND DISSPLA NEEDS ADDITIONAL WORKSPACE 
BUR SMOOTHING. SAME COMMENTS CORRESPONDS 10 THE 
peatll BCOMON" STATEMENT. 


COMMON WORK(5000) 
Se) PARAMETERS FOR AXES: 


DEFINE X 


AXIS PARAMETERS 


XORIG=0.0541 


XSTP= 


DoS 


XMAX=2.9541 


DEFINE Y AXIS PARAMETERS CCORRESPOND TO 2 AXIS) 
NORIG=0\.115 
Mowe =0 256 
YMAX=3.015 
Seer COMPRS 10 USE THE DISSPOP OPTION. 
CALL COMPRS 
GALE SETCLRC* CYAN") 
Semernce AND PLOT SIZES, SET UP AXES FOR PLOT: 
GAPE TPAGECS.. 511.0) 
CALL BCOMON(5000) 
CALL AREA2D(6.0,8.0) 
LABEL AXES: 
CALL XNAMEC'X-AXIS VARIATION CMETERS)$',100) 
CALL YNAMEC'Z-AXIS VARIATION (METERS)$*',100) 
CALL GRAFCXORIG, XSTP,XMAX, YORIG, YSTP, YMAX) 
CALL FRAME 
fenee OF THE PLOT: 
Cael HEADING "CONTOUR E-FIELD (DB REF [0 1V7MJis*,100,2.,3) 
CAEL HEADINC’MONOPOLE 6CM AT CORNER OF WIRE GRID BOX$',100,1.,3) 
CALL HEADINC *FREQ=1GHZ-(5.3 CM ON DIAGONAL)$',100,1.,3) 


MAKE CONTOURS AND DRAW: 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


END 


SEECLERC' RED") 

CONMIN(3.0) 

CONANG(60.) 
CONLINCO, 'MYCON', 'LABELS',2,10) 
CONMAKCA,NX,NY,FINC) 

CONTUR(1, 'LABELS', 'DRAW') 


ENDPL(O) 


Cc 
CHK HHH KH HK HH HH HK HK HK HK HH HK HH HK HK HH HHH HK HHH HK HHH HH HH KKK HH HK HK KKH KKK KKK KK KKK 


C 
C THIS ROUTINE MAKES NEGATIVE CONTOURS DASHED , 


SUBROUTINE MYCONCRARAY, IARAY) 
DIMENSION RARAY(2),IARAY(1) 


ZERO LINES HEAVIER 


17a, 


CON00010 
CONO0020 
CONO0030 
CON00040 
CONO00050 
CoOn0d0060 
CON00070 
CON00030 
CON00090 
CON00100 
CONO0110 
CON00120 
CON00130 
CON00140 
CON00150 
CON00160 
CON00170 
CON001380 
CON00190 
CON00200 
CON00210 
CONOO220 
CON00230 
CONO00240 
CON00250 
CON00260 
CON00270 
ConNddg2z30 
CON00290 
CONO00300 
CON00310 
CON00320 
CON00350 
CON00340 
CON00350 
CONO0360 
CON00370 
CONO0330 
CON00390 
CON00400 
CON00410 
CON00420 
CON00430 
CON00440 
CON00450 
CONO0460 
CON00470 
CON004380 
CON00490 
CON00500 
CON00510 
CON00520 
CON00530 
CON00540 
CON00550 
CON00560 
CON00570 
CON00580 
CON00590 
CON00600 
CON00610 
CON00620 
CON00630 
CON00640 
CON00650 
CON00660 
CON00670 
CON00630 
CON00690 
CoON00700 
ConN00710 
CON00720 


Pree 


CONTOR FORTRAN Bl 


C WHILE THE POSITIVE CONTOURS ARE PRESENTEDSBY SOLID Lies: 
c 


10 


CALL RESETC’ DASH") 

IF CRARAYCI) GE. G2) GO" TO 20 
CALL DASH 
RARAYC2) = 
IARAYC1) = 
IF CRARAYC 
RETURN 

END 


1 
1 
1) .EQ. OG.) TARAYCI) = 2 


CONTOR FORTRAN: Subroutine Which is Used to 
Generate the Contour Plots. 
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CON007 30 
CON00740 
CON00750 
CON00760 
CON00770 
CON00780 
CON00790 
CON00800 
CON00810 
CON00820 


Elle: CONS FORTRAN Bl 


© thIS PROGRAM READS THE VALUES OF THE PHASE OF EX OR EZ COMPONENT 
ESEROM A NEC OUTPUT FILE Glee OUTPUT IN THIS CASE HAS TO 
eect PRODUCED BY A PU 2r2s¢e) ailb Ne CARDS .THE QUTPUT 
Peeeevits CASc 15 LIKE THEeOile Ww EXAMPLE 2 DATA FILE) 
C AND USING THE SUBROUTINE CONTOR FORTRAN PLOTS THE 
: VARIATION OF THE PHASE USING DISSPLA. 

REAL*8& RO 
PReOEPINE THE ARRAY ACNXsNY) WHICH IN THIS CASE IS COMPOSED 
C OF 100X100=10000 ELEMENTS. NX=NUMBER OF COLUMNS AND 
Cehy-NUMBER OF ROWS OF THE ARRAY. 

DIMENSION AC100,100) 

EMIN=1E20 

EMAX=-1E20 


DO 50 I=1,100 
DO 40 J=1,100 
READ €12,20)R0 
BeUSE THIS FORMAT HHEN THE PHASE OF X COMPONENT 
3 Uiernie ELECTRIC FIELD IS 70 BE FLOTPED 
¢ 


20 FORMAT €25X,D10.3) 
feeoc THIS FORMAT WHEN THE PHASE OF Z COMPONENT 
Secnetne ELECTRIC IS TO BE PLOTTED 
C 


C20 FORMAT (€65X,D10.3) 

ACJ,I)=RO 
Pern tHIS STEP THE PROGRAM PREDICTS THE MINIMUM 
3 AND MAXIMUM OF THE PHASE VALUES. 


EMIN=MINCEMIN,ACJ,I)) 
EMAX=MAXCEMAX,ACJ,19) 


40 CONTINUE 
50 CONTINUE 
NxX=100 
NY=100 
BeUEFINE THE NUMBER OF LEVELS BETWEEN CONTOURS 
NLEVEL=10 
EOIN THIS STEP THE PROGRAM CALCULATES THE CONTOURING 
4 INTERVAL 
PINC=( EMAX=EMINIZCNLEVEL=1 ) 
C THE VALUES OF EMIN,EMAX,FINC ARE DISPLAYED ON THE 
C SCREEN 
Cc 
WRITE €6,%) EMIN,EMAX, FINC 
Caine CALLING STATMENT GF THE SUBROUTINE CONTOR FORTRAN 
CALL CONTOR CA,NX,NY,FINC) 
STOP 
END 


6. PROGRAM II: CON3 FORTRAN. 
for Near Electric Field. 


Contour Plots of Phase 


Pele: CONS EXEC Bl 


FILEDEF 12 DISK FILENAME PLOTDATA A 
EXEC DBISSPLA CONS 


fen CONS EXEC: 
of DISSPLA. 


Fortran File for Automatic Execution 
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Cornoool1o 
CONO0020 
Conoooso 
CON00040 
CONO0050 
CoNoooéC 
CONO0070 
CONd0020 
CoNoo090 
Cornoo100 
Conogo110 
Conoo120 
CONoo1S0 
CON00140 
CONO00150 
CONO0160 
CONO00170 
CON00180 
CON00190 
ConNo0o0z00 
CONCO2Z10 
Conoo220 
CON0OQ230 
CON00240 
CONO0250 
COH00260 
CON00270 
CON002380 
CoNd0290 
CoNoo0300 
CONO0310 
Conoog320 
CON00330 
CoNo0340 
CONooz50 
CON00360 
CON00370 
CON00380 
CON00390 
CON00400 
CON00410 
CON00420 
CON00430 
CON004490 
CON00450 
CON00460 
CON00470 
CON00480 
CON00490 


FILE: CONZ FORTRAN Bl 


THIS PROGRAM READS THE VALUES OF THE MAGNITUDE EREAR 0G e- TEED 
FROM NEC OUTPUT FILE AND GENERATES A PARTICULAR PLOT WHICH HAS 
THE 3-D SURFACE AT THE TOP WHILE AT THE BOTTOM THE SCONmOUR 

PLOT IS PRESENTED. THE PROGRAM USES THE SUBROUTINES Di HGg ZU scr 
DISSPLA GRAPHICS PACKAGE WHICH IS ATTACHED HERE AT THE END. 


ANQNANANO 


REAL*®8& RO 
DIMENSION AC30,30) 
CALL COMPRS 
EMIN=1E20 
EMAX=-1E20 
DO 50 I[=1,30 
DO 40 J=1,30 
READ (12,20)R0 
20 FORMAT (€17X,D10.3) 
ACJ,1)=20*LOG10CRO) 
C IN THIS STEP THE PROGRAM PREDICTS THE MAXIMUM AND MINIMUM VALUES 
C OF THE MAGNITUDE OF THE ELECTRIC FIELD: 


Cc 
EMIN=MINCEMIN,ACJ,1)) 
EMAX=MAXCEMAX,ACJ,1)) 
40 CONTINUE 
50 CONTINUE 


NX=30 

NY=30 
C DEFINE THE NUMBER OF LEVELS CINPUT FROM KEYBOARD) 
Cc 


WRITE (6.70) ENTERSEPE VEER a 
READ C5y x2) NLEVER 
NLEVEL=20 
IN THIS STEP THE PROGRAM CALCULATES THE CONTOURING INTERVAL. 
FINC=C EMAX-EMIN)/ CNLEVEL=-1) 
THE VALUES OF EMIN,EMAX,FINC ARE DISSPLAYED ON THE SCREEN. 
WRITE €6,%) EMIN, EMAX,FINC 
IN THIS STEP THE USER HAS TO INPUT BY THE KEYBOARD THE COORDINATES 
OF THE VIEW FONT. 


WRITE €6,70) * ENTER VA5VZ,VE.> * 
READ (€5,%) VX,VZ,VE 
THIS IS THE CALLING STATEMENT OF THE SUBROUTINE "DMH0O2™ 


CALL DMHOO2 CA,NX,NY,FINC, EMIN-30, EMAX, NLEVEL,VX,VZ,VE) 
CALL DONEPL 

70 FORMAT CA) 
21 OP 


END 


ga0N90N OQ 9A 


an 


CHK HK KKK HK HK HK HH HHH HK HHH HH HHH HH HHH HH HH HH HH HHH HH HHH HH HHH HHH HH HH HK HK KK KH KKK KK KKK 


SUBROUTINE DMHOO2CZDAT,NX,NY,FINC, EMIN, EMAX, NLEVEL,VX,VZ, VE) 
DIMENSION ZDATC30,30) 
COMMON WORK(C5000) 
CALL RESET (SHALL) 
CALL (PAGE Cli 55:5).5) 
C SET AXIS PARAMETERS AND ALPHABETS 
CALL INTAXS 
CALL ZAXANG (90.) 
CALL SIMPLX 
CALL MXLALE CSHEZCST, 22°) 
CALL MX2ALF CSHSTAND,'%') 
C DEFINE 3-D WORK AREA AND AXIS 
CALL AREAZD (9.5,6.5) 
C DEFINE THE LABELSOF THE Eom 
CALL HEADIN ('%3-D ELECTRIC FIELD (DB REF TO 1V/M)$',36,2.0,3) 
hee HEADINC'%MONOPOLE 6CM AT CORNER OF WIRE GRID BOX$',40,2.0, 


CALL HEADINC'%(5.3 CM ON DIAGONAL FROM CENTER)$',33,1.0,3) 
CALL RESET C6HHEIGHT) 
C DEFINE THE LABELS FOR THE THREE AXES 
CALL X3NAME C'%X VARIATION (CMETERS)$',21) 
CALL Y3NAME (€'%Z VARIATION (METERS)$',21) 
CALL ZSNAME C"ZELECTRIC FIEUDSIC TS) 
CALE VOLMSD (10.,107,56.) 
C DEFINE 3-D VIEW POINT 
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CON00010 
CON00020 
CON00030 
CON00040 
CON00050 
COHO0060 
CON00070 
CON00080 
CON00090 
CON00100 
CON00110 
COHO00120 
CON00130 
CON00140 
CON00150 
CON00160 
CON00170 
CON001380 
CON00190 
CON00200 
CONO00210 
CON00220 
CON00230 
CON00240 
CON00250 
CONO0260 
COHN00270 
COH002380 
CON00290 
CON00300 
CON00310 
CON00320 
CON00330 
CON00340 
CON00350 
CON00360 
CON00370 
CON00330 
CON00390 
COHN00400 
CON00410 
CON00420 
CON00430 
CON00440 
CON00450 
CON00460 
CON00470 
CON004380 
CON00490 
CON00500 
CON00510 
CON00520 
CON00530 
CON00540 
CON00550 
CON00560 
CON00570 
CON00580 
CON00590 
CON00600 
CON00610 
CON00620 
CON00630 
CON00640 
CON00650 
CON00660 
CON00670 
CON00630 
CON00690 
CON00700 
CON00710 
CON00720 


FILE: CON2 FORTRAN Bl 


CALL VIEW (VX,VZ,VE) 


F CALL GRAF3DCXORIG, XSTEP, XMAX, YORIG, YSTEP, YMAX, EMIN, FINC, EMAX) 
“ar la C020561,.5252.9591,0.115,.53,3.015,EMIN, FINC, 
xX 
C DEFINE THE SURFACE VIA A MATRIX 
CALL BLSUR 
CALL SURMAT (ZDAT,2,30,2,30,0) 
CALL DASH 
XDEL=0. 


DOV 2si T1535 
CALL BSHIFT (XDEL, 0%) 
C DRAW BEDPOST EFFECT 
CALL RESETC'BASALF®) 
CAEL HEIGHT. 2) 
CALL COMPLX 
Se DEFINE A MESSAGE ON THE PLOT IF IT IS DESIRED 
c CALL MESSAGC'SUBROUTINE DMHO02$!',100,3.5,0.) 
CALL RESETC*HEIGHT ') 
CALE RESETC *COMPLK™) 
DEFINE THE 4 DASHED LINES WHICH CONNECT THE 4 CORNERS OF THE CONTOUR 
GRAPH WITH THE 3-D GRAPH.THESE 4 DASHED LINES ARE AN OPTIONAL FEATURE. 


CALL RLVECS (0.0561,0. 115, ZDAIC1,1),0.0541,0-.115, EMIN,0) 
CALE RLVECS C2.9561,0.1715,2ZDA1(50,1),2.9591,0-115,EMIN, 0) 
CALL RLVECS (€2.9541,3.015,ZDAT(30,30),2.9541,3.015,EMIN,0) 
CALL REVECS (0.056153. 015;,ZDAT( 1,30), 0. 0541,3.015, EMIN, 0) 
XDEL=XDEL+.01 
eo.) CONTINUE 

CALE RESET (SnD0T) 
CALL RESET (6HBSHIFT) 

SeeenteR GREFITI LOOP AND DEFINE 2=D PLOT 
CALL GREITI CO..50).50.051.5,60.50.050771 -,0..) 
CALL INTAXS 
CALL AREAZD (10-,10.) 

CALL GRAF (CXORIG,XSTEP,XMAX,YORIG, YSTEP, YMAX) 


CALL GRAF €0.0541, .5652.9541,0. 115 -58,3.015) 
C CONTOURING IS OPTIONAL FEATURE AND CORRESPONDS TO THE BOTTOM PLOT 
C SET CONTUR PARAMETERS 
CALL BCOMON (5000) 
CALL CONMAK C(ZDAT,30,30,FINC) 
c CALL CONLIN (€4,5HSOLID,6HLABELS,2,10 
CALE CONLIN (0, "MYCON", "LABELS", 2, 10 
CALL CONANG (90.) 
CALL FRAME 
SCALU HEIGHT C.3 
Cc CALL CONTUR (5 
CALL CONTUR (1 
END GRFITI LOOP 
CALL END3GR(0) 
CAEL ENDPL (0) 
RETURN 
END 


C 
CH HH KK HHH HH HH HH HK HH HHH HH HHH KK HH KH HHH HK HHH HK HH HK HHH HH HH HK HH KH HK KKK KKK KKK 
Cc 


aAN0 


) 
) 


) 
6HLABELS, 4HDRAW) 
"LABELS", 'DRAW') 


, 
, 


oO 


SUBROUTINE MYCONCRARAY, IARAY) 
DIMENSION RARAY(2), IARAY(1) 


Cc 
C THIS ROUTINE MAKES NEGATIVE CONTOURS DASHED AND THE ZERO LINE HEAVIER. 
GeuHE POSITIVE CONTOURS ARE PRESENTED BY SOLID LINES 
CALL RESETC'DASH') 
IF CRARAY(1)- .GE. 0.) GO TO 10 
CALL DASH 
10 RARAY(2) = 1. 
IARAY(1) = 1 
IF CRARAYC1) 
RETURN 
END 


1) .EQ. 0.0 TARA Ge =—2 


8. PROGRAM III: CON2 FORTRAN. 3-Dimensional Plots of 
Magnitude for Near Electric Field. 
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CON00730 
CON00740 
CON00750 
CON00760 
CON00770 
CON007380 
CON00790 
CON00800 
CON00810 
CON00820 
CON00830 
CON00840 
CON00850 
CON00860 
CON00870 
CONO00S80 
CON00890 
CON00900 
CON00910 
CON00920 
CON00930 
CON00940 
COHN00950 
CON00960 
CON00970 
CON00980 
COnN00990 
CONO1000 
CON01010 
CON01020 
CONO01030 
CON01040 
CON01050 
CON01060 
CON01070 


Oo 

oO 

mae 

oO 

ko 

— 
WOONAtINEAWGAYH Ow O 
e0oqooo0o oo 0oo0c a0 0 


CON01200 
CON01210 
CON01220 
CON01230 
CON01240 
CON01250 
CON01260 
CON01270 
CON01230 
CON01290 
CON01300 
CON01310 
CON01320 
CON01330 
CON01340 
CON01350 
CON01360 
CON01370 
CON013380 
CON01390 
CON01400 
CON01410 
CON01420 


Febes (cone EXEC Bl 


FILESES l2 DISK FILENAME Poon gatas 
EXES DISSPLA CONZ2 
9. CON2 EXEC: Fortran File forsAutonat 7. 
Execution of DISSPLA. 
PART B: NEC INPUT DATA FOR NEAR FIELDS. 
SURFACE PATCH MODELING. 
FILE: DIPS DATA Bl 
CM DIPOLE L’2 ON Z AXIS.FREQUENCY 1 GHZ 
= IN FREE SPACE. NEAR E=rieLD CS0xKS0M rors? 
GH 1,5,0,0,-0.075,0,0,0.075,050016 
Ge 
FR 0,0,0,0,1000 
PMO lee ld 
Pie ceo 
NE 0,30,1,50,0.0l66,0,0.115,0.1 000 
XQ 0 
Er 
1. Data File Dipole A/2 in Free Space. Near 
Electric Field (for Contour Plot of Magnitnd ge 
FILE: DIPHS4 DATA Bl 
C4 DIPOLE L7Z2 ON Z AXIS FREQUENCY Se lGageuierReaese eer. 
CE NEAR E-FIELD.100X100 POINTS. 
ee 1,5,0,0,-0.075,0, 070607 5,0..0016 
ic 
FR 0,0,0,0,1000 
MoO pls Oued 
are ar sarah 
NE 0,100,1,100,0.0166,70,0.115,70-029,0,0 029 
XQ 0 
a EN 
2. Data File Dipole A/2 in Free Space. Near 


Electric Field (for Contour Plot ei Phase): 


TS2 


PLE: MONA D1 DATA Bl 


C— OVER PERFECT GROUND. NEAR E-FIELD (30X30 POINTS). 

Gil Ss O00, OO), O06) 0 001.6 

GE l 

Gtr i 

FR 0,0,0-0, 1000 

 Wotlotl oe niyo 

iy Bawaba 

HE Boo os ee 

aA 

EN 

3. Data File 6cm Monopole Over Ground Plane. 

Near Electric Field (for Contour Plot of 
Magnitude). 

FILE: MONPHl DATA Bl 
CM MONOPOLE 6CM ON Z AXIS.FREQUENCY 1 GHZ 

CE OVER PERFECT GROUND. NEAR E-FIELD.C100X100 POINTS). 
Gi 2-5,0,0,0,050,0.06,0.0016 

GE l 
GN 1 

Ie 40s Eel alto 6) 

Pee, OS 0 

PL 2,2,4,1 

HE 0,100,12,100,0.0166,0,0-175,.0 2029, 0,0.029 

XQ 0 

EN 

4. Data File 6cm Monopole Over Ground Plane. 


MONOPOLE 6CM ON Z AXIS.FREQUENCY 1 GHZ 


Near Blectric Field (for Contour Plot of 
Phase). 
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BLES DATA Bl 


CH SPEOXK AP PA CHES ON Or SUREACE 

CE FREQUEHCY 21.00 GHZ 

SHS550 soo Uo 0 on Ue 

SM35 Sp 0.55 -0eo. U- Ue, 0s 0 

SO Os SO 

GX 0,110 

SM7 5 7770.55-0. 5515025. 0 oe 

3€0,0,0 95,0" 552 

GS 0,0;0.1 

GE ol 

FR 0,0,0,0,1000 

GN ol 

VIG 

HX 

CM 6 CM MONOPOLE AT CENTER GF PATCH BOX 230X300 POrNns 
Ee NEAR E-FIELD ON X-AXIS AT 0.0166,0,0.115(CONTOURED FIELD) 
Beye nec at ee ame 


1 
30,0.0166, 0,0 .115,021,0, 021 


> 
Q2Onoro 


5. Data File Gem Monopole at Center of Suriace 
Patch Box. Near Electric Field (for Contour 
Plot of Magnitude). 


FILE? NEAL DATA Bl 


CM SPBOX 7X7 PATCHES ON TOP SURFACE 

CE FREQUENCY 1 00yGizZ 

S13 5570.57 70.5;0;0.5,-0., 571 

SICUs O05 0a sort 

SiS 5 Ol oe Or OO >a Cisd0) 

SC0me, 025,052 

GX 0,110 

SiW7 77 3-0. o> 02551, ee Oy Oe oe 

SCO .0, 05,0. 5,1 

GS.0,0 071 

GE 2 

FR 0,0,0,0,1000 

GN 1 

VIG 

NX 

CM 6 CM MONOPOLE AT CENTER OF PATCH. 20x 
CE NEAR E-FIELD ON X-AXIS AT 0.0166,0,0.115.POINTS 100X100 


GF 
GH1,5 
GS 


? 


ie aie ca 


? 


21,0 


0 
1 1 
2 

100,1,100,0.0166,0,0.115,0.029,0,0.029 


0 
71,0 
pias: 
00,1 


6. Data File 6cm Monopole at Center of Surface Patch 
Box. Near Electric Field (for Convour ore 
Phase). 184 


Rabe: NED! DATA Bl 


CM 11X11 PATCHES ON TOP SURFACE (6CM AT EDGE 3.63CM FROM CENTER) 
CE FREQUENCY 1.00 GHZ 

Sains 0 Se 6 os Or Oe oe culos 

SCO; 0, Oy - 0). 5,1 

SM3,5,0. CS SU) Sin U 10 

$C0,0,0.5,0,1 

GX 0,110 

Sills O0. 5,-0.4 55,0. oy Uno 1 

Se0,0,0.5,0.5,1 

Gs 0,0,0.1 


FR 0,0,0,0,1000 


Q 
= EH 
ror 


CM 6 CM MONOPOLE AT EDGE 3.63 CM FROM CENTER (11X11 PATCHES ON TOP) 
CE NEAR E-FIELD ON X-AXIS (CONTOURS 30X30 POINTS) 


GIN 1,5,70.363636,0,1,-0.363636,0,1.6,0.016 


l 
30,—-0,0197656,0.0.115,0,1,0,0.1 


me 
or 
OONOrFOS 


7. Data File 6cm Monopole at Edge of Surface Patch 
BOx. Near Electrre Field "tor scontour Plot of 
Magnitude). 


Bie wer H) DATA Bl 


GM Lixll PATCHES ON TOP SURFACE C6Gi AT EDGE 3.63CHM FROM CENTER) 
SCE FREQUENCY 1.00 CHa 

Sool, 0s 55 0, 0, 5,20. So 

SCO, 0, 0,-0.5, 1 

Silas, 05s -0.5,0,0.5,0,0 

560,0,0.5,0,1 

GX 0,110 

Si eleeleles Ol 55 — Ole Dip lip One Dip Oe oe 

S605 050,5,0.5,1 

GS 0,0,0.1 


a 
FR 0,0,0,0,1000 
i 


Ciiltowch MONOPOLE AT EDGE 3.635 CM EROS GENTERSCEIXII PATCHES ON TOP) 
Ge NEAR E-FIELD ON X=AXIS CPHASE 100 PUINTS-ONE NE CARD) 


Gh 1,5,-0.363636,0,1,70.363636,0,1.6,0.016 
GS 105.0,0.1 


0 
00,-0.0197656,0,.0. 115, 070297510, .0-7 029 


8. Data File 6cm Monopole at Edge of Surface Patch 
Box. Near HlectricelteldsctonscContour Plot of 
Phase). Mees, 


FILE: NeCL DATA Bl 


CM 11X11 PATCHES ON TOP SURFACE (6CM MONOPOEE Ad COmiER S714 7¢7) 
CON THE DIAGONAL FROM CENTER 

CE FREQUENCY 1 .00meGHZ 

Spe) 5 Fn Ula sin Oy Wo Sin ll 6 Bip Jk 

SG, OG.0 — Oe ok 

$13,353,005) -0 Su0n0ns 020 

5C0,50,,075, 071 

GX OF710 

SMI 3 170255 =0.55., 0. oo— ce oe. 

SCO, 05.0. 5.0.55 1 

G5 7050500 

GEA 

Fr 0, 0,9,0,1000 

GN l 

WG 

NX 

CM 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
Cii TOP SURFACE L1IX11 -PATGHES 

Cit FREQUENCY Si 00 3GHz 

CE NEAR E-FIELD ON X=AXIS2CCONTOURKS 30X30 POINTS) 

GF 
GH1,5,0'.363636,0.363636,1,0.363656,0.365656,1 6, 0.046 


pls 
50,70..0529636,7 0036565670 -115,051 700 


9. Data File 6cm Monopole at Corner of Surface Patch 
Box. Near Electric Field (for Contour Plot ws 
Magnitude). 


FILEs NGOPHI DATA Ba 


Ci 11X11 PATCHES ON TOP SURFACE (6CM MONOPOLE AT GERHER 5014.cup 
CM ON THE DIAGONAL FROM CENTER 
CE SPREQUENGY 1700) GHZ 
HSp 3 Olea. — Und, 0p leo — 0). 5 
SG0, 0005, 1 
Sitoy O5.Un os Us 5p UZ Uo, Ua 
SCO, 0,055,071 
GX 0,110 
SM ijl = 0's 0. Oo, Oo See 
5C0;, 05.0.5, 07.5, 1 
GS 0, 0,0. 
Geri 
FR 0,0,0,0,1000 
al 


Ci’ 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
CM TOP SURFACE 11X11 oPAgCHES 

CM FREQUENCY 1.00 GHZ 

CE NEAR E-FIELD ON X-AXIS.CPHASE 100 POINTS-ONE NE CARD) 


GF 
oy »5,0.363636,0.363656,1,0.363636,0.363636,1.6,0.016 


S05 05.0.1 
GE l 
x Opt UG ru 
PL 2,2,4,1 
ie poe eee 2 0.0529656,0.0565656,0.115,0-029,,050. 0c 


10. Data File 6cm Monopole at Corner of Surface Patch 
Box. Near Electric Field (for Contour Plot of 
Phase). 186 


PART C: RADIATION PATTERNS FOR MONOPOLE 
ON THE SURFACE PATCH BOX. 


FILE: SPBOXV1 DATA Bl 


CM SPBOX 7X7 PATCHES ON TOP SURFACE 
CE FREQUENCY 1.00 GHZ CVERTICAL PATTERN) 
Ss 0.57-0.5)0,0-5,-0.5,1 
Seu, 070, 0.5, 2 
Ss opoe 0. 55-0. 5,050.55 070 
S60 ,0,0.5, 0, 1 
GX 0,110 
ST lee OMS pit Ves 30 (bps) 20) ores 
Si OU 5 OS ak 
657707070. 2 
Ge 1 
FR 0,0,0,0,1000 
ik 


Ci4 6 CM MONOPOLE AT CENTER OF SURFACE PATCH FREQ. 1.90 GHZ 


GM FOP SURFACE 7X? PATCHES 
E RADIATION PATTERNS CVERTICAL) 


GF 
GW1,5,0,0,1,0,0,1.6,0.016 


G570;,0,0.2 

CE 

Exd0> lod, 0,159 

PL 3,1,0,4% 

RP OF; lal, 1,1002,-9950,.1,0 
XQ 0 

EN 


1. Data File 6cm Monopole at Center of Surface Patch 
Box. Vertical Pattern. 


FILE: SPBOXH1 DATA Bl 


Chase pOx fr? PATCHES ONSTOP SURFACE 

CE FREQUENCY 1.00 GHZ CHORIZONTAL PATTERN) 
Sms n Sip Ueno as Wh Sp S05 baal 

Se Oils 0,-0.5, 1 

Sito Ul, Ul. o7 0s Os oe U5 0 

Se0, 0,0. 5,071 

Gre 07 2.0 

Srl ot po Oe Ol ool 0. Oe 

SC0,0,0,.0> 0.251 

G520;5,0;0.2 

GEat 

FR 0,0,0,0,1000 

GN 1 

WG 

HX 

CM 6 CM MONOPOLE AT CENTER OF SURFACE PATCH FREQ. 1.00 GHZ 
C4 TOP SURFACE 7X7 PATCHES 

CE RADIATION PATTERNS CHORIZONTAL) 


Gr 
GW1,5,0,0,1,0,0,1.6,0.016 


, 


al 
»1002,90,0,0,1 


0,0. 
2 a 
Uy 4. 
pl sol 


2. Data File 6cm Monopole at Center of Surface Patch 


Box. Horizontal Pattern. 
187 


FILE: SPEDGEY) Dane Bl 


CM SPEDGE1] 11X11 PATCHES ON TOP SURFACE (6CM AT EDGE 3.63CM FROM CENTER) 
CE FREQUENCY 1.00 GHZ (VERTICAL PATTI ERH-xX fais CUD 

SUA on Uo Sin —tle Sn nl) 6 5p SO, Soil 

SCU, 0,0. =07 5c. 

SiWei5 Sin Uo Bn OO Sn Wn San 

$C0,070.5,050 


GX 0,110 

SM, ll 02s s oe lo ore 
SCO, U5 Uo oote on 

GS. 0,0,0.1 

GE 

FR 0,0,0,0,1000 

GN 1 

HG 


Hx 

Ci 6 CM MONOPOLE AT EDGE 3.63CM FROM CENTER FREQUENCY 1700 GHz 
Ci! TOP SURFACE lik? PATCHES 

CE RADIATION PATTERNS CVERTICAL-X AXIS CUT) 


GF 

GW1,5,-0.363636,0,1,-0). 36565650, 1.6..0).006 
G5 0;0 7072 

GE 

EX 0,316 071,0 

PL 3515054 

RP 0,161,1, 1002, - 70,0, 8 

XQ 0 

EN 


3. Data File 6cm Monopole at Edge of Surface Patch Box. 
Vertical Pattern (x-Axis Cut, 


FILE: SPEDGEVZ DATA Bl 


CM SPEDGE] 11X11 PATCHES ON TOP SURFACE <6CM AT EDGE 3°63CiH FROM Cee 
CE FREQUENCY (2. GHZ Se ee PATTERH-Y AAIS CUdD 

SC0, 0,0) -0ns sl 

SMS 570.5. 20 .5,50,0.5,0,70 

5C0,0,0.5, 07.1 
GX. 0,110 

SHI 1s 30 5550.55.150.5.50 55.1 
SCO,0,0. S50c Sel 

Gs 0,0,0. 
GE I 

FR 0,0,0 
GN l 
WG 

NX 

CM 6 CM MONOPOLE AT EDGE 3.63CM FROM CENTER FREQUENCY 1.00 GHZ 
Ci -TOP SURFACE TIX] DP PATCHES 

ce RADIATION PATTERNS CVERTICAL-Y AXIS CUT) 
GW1,5,=0.363636,0,1,-0. 363565650, .c70 016 

GS 0,0,0 


1 
»0,1000 


= 
— 


002,-90,90,1,0 


4. Data File 6cm Monopole at Edge of Surface Patch Box. 
Vertical Pattern (y-Axis Cut). 
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Bl 


Ti] 
=] 
- 
ma 
“” 
U 
nm 
u) 
@ 
n 
| & 
hr 
oO 
J> 
Ho 
a 


Dap se oe oe ero resUcrEG. S6C Al —EXGE S,65CM FROM CENTER) 
00 GHZ CHORIZONTAL PATTERN) 

pA). Be 18) ,2}p dl 
OF 


=) EG) 
1 


intin tinny 
iG 


~~ ©» 


HX 
Gio CM MOHOPOLE AT EXGETS.65Cm FROM CENTER FREQUE:ICY 1.00 GHZ 


BaemOr SURFACE 11X11 PAREHES 
CE RADIATION PATTERNS CHORIZONTAL ) 


Gr 

Gil, 5,-0.565626,0,1,=0.565656,0,1.650.016 
GSy0),o. OF 2 

Ge 

Pi Ls Ly Ug >, O 

Be ores 0.4 

Re O,1,361,1002,90,0,0,1 

xQ 0 

EN 


o. Data File 6cm Monopole at Edge of Surface Patch Box. 
Horizontal Pattern. 


= 
yY> 
Y- 
od 
t 


meee: SPCY 1 


C4 SPCNTR 11X11 PATCHES ON TOP SURFACE (6CM MONOPOLE AT CORNER 5.146 CM 
CaO The DIAGONAL FROM CENTER 

Seen elsUciCyY 1.00 GHZ. VERTICAL PAITERN=Cx-AXIS CUT) 
Ss, Sy 0. Sp SU Se Ula Se Sle al 

Sco 0gl,—0.5, 1 

Ties, 0.5570.5,0,0.5,0;,0 

seo, 0,0.,5, 0, 1 

GX 0,110 

Sil, el ls 5, =O, Syl, Oe =O). 5, 1- 

Sseo,0;,0,5,0.5,1 

cS 0,050.1 

GE] 

FR 0,0,0,0,1000 

GN 1 

e 

IX 

Cit 6 Ci MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
Gaelor gsUREACE LIXll PATCHES 

Gy FREMUEHCY 1.00 GHZ 

Ce RaUte LON PATTERNS CVERTICAL= X AXIS CUT) 


GF 
GW1,5,0.363636,0.3636356,1,0.3636356,0.3635636,1.6,0.016 

Gs 0,0, 0).1 

Ge t 

PO, 2,051, 

PE s,1,0,49% 

Ree oe 002, -90,0,1,0 

XQ 0 

EN 

6. Data File 6cm Monopole at Corner of Surface Patch Box 


Vertical Pattern €(x-Axis Cut). 


Leg 


FILES SEGV2Z DATA Bl 


CM SPCNTR 11X11 PATCHES OW TOP SURFACE C6CM MONOPOLE Al CORNER 5. 1¢ecy 
Cid ON THE DIAGONAL FROM CENTER) 

CE FREQUENCY 1.090 GHZ. VERTICAL PATTERN <95) DECKEES Culp 
Si35 550). 57-0855 OOo Oreo nl 

SC05 0,0, 0-550 

Siise 0.0, 0.57070 so Ono 

5C0;0;0.5, 01 

GX 0,110 

Sie le — Ooi 0S Oo Oro 

SC0,0, 0.57 0051 

G50,0;, 0.1 

Gel 

FR 0,0,0,0,1000 

GH 1 

WiG 

NX 

CM 6 CM MONOPOLE AT CORNER 5.19 CM ON THE DIAGONAL EROM@GENTER 

CM TOP SURFACE 11X17) PATCHES 

CM FREQUENCY 1.00 GHZ 

CE RADIATION PATTERNS. VERTICAL €45 DEGREES CUT) 


GF 
Ghil 55> 0:.363636,0,35365656,1,0., 526565670 7505050. oer une 
iI 


GS 0,0,0. 

Geni 

x Fl 0; 10 

PL 3,1,0,4% 

RP 0,181,1,1002,-90,45,1,0 
XQ 0 

EN 


7. Data File 6cm Monopole at Corner of Surface Patch 
Box... Vertical Patterm (tor stut). 


FILES sea DATA Bl 


CM SPCNTR 11X11 PATCHES ON TOP SURFACE €6CM MONOPOLE AT CORNER 5.14 CM 
CM ON THE DIAGONAL FROM CENTER 

CE OPREGUENCY 200 SGHZ. HORIZONTAL PATTERN 
SiSeoe leo reo. 0. Un o> Oleos 

$C0,0,0>-0.5,1 

Sitse 550 54 Ooo, U0. 5,050 

500,050 5,0. 1 

GX 0,110 

SM1L1,11,—-0.5,-0.5,1,0.5,-90.5,1 

SC0707.07550'. 5,1 

GS 0,0,0.1 

Gea 

FR 0,0,0,0,1000 

GN 1 

N 

Ci4 6 CM MONOPOLE AT CORNER 5.14 CM ON THE DIAGONAL FROM CENTER 
Ci TOP SURFACE TIX11 PATEHES 

CM FREQUENCY 1.00 GHZ 

CE RADIATION PATTERNS CHORIZONTAL) 

GF 
GW1,5,0.363636,0.363636,1,0.363636,0.363636,1.6,0.016 


002,90,0,0,1 


8. Data File 6cm Monopole at Corner of Surface Patch 
Box. Horizontal Pattern: 


coe 


APPENDIX C 
NEC INPUT DATA FILES 


PART A: NEC INPUT DATA FILES FOR WIRE GRID 


MODELING. STRUCTURE GEOMETRY 


PLLEs WG DATA Bl 
CM WIRE GRID 10 CM. LONG WITH 6CM MONOGPOLE AT CENTER 
CHOCELLS pilesexk «leo OROuHD 
CE DATA FILE FOX PLOVIING UsING PEOTRDGLP EXEC 
Gh LOO, Gp aay eee le pal erik 
Gif 1.7,0;0,0;0, 0395125, 100.100 
GH ZNO A 5 6 Sno Bodin co Spa Ss Qn o Wk 
Gri 2 An Us On lWplin o LAS, 0, ZUG. Zeal 
GM 55 5o 05s 0, >. eoro, UrcuU. col 
GHW Ss Min 5 SMG Sin ln oS Wo he a Ot 
Gi4 CAT nn Op On pW po, Ley SUG) 6 SOS 
Gil G00, 4505. . oF cele te. 51 o> ue 
Gif 5,2,0,0,0,0,.125,0,900.404 
GH S006) ..57 5.—. 56h, 57 Se Uo, 20d 
Gi G52,070,0,—-.les.0, 07 500.505 
Gx sun ee 
GH 99,6... 5,0,0,.o70n1, won 
Gil f> so Ores 0 Oren. eae 
Gri COs ~aoe lo Us poe oe ed 
Ghi Seon coe le ly so, Deda. Ow 
Gh 1, 5,0.005,0,1,0.005,0, 13654076 
4 0,35,0,0,10,0,0,0,001.001 
GE 
Eit 
1. Data File Geometry of 6cm Monopole at Center 
Ot WirevGrid box. 
FELE Nelo DATA Bl 
Ci WIRE GRID 10 CM. LONG WITH 6CM MONOPOLE AT EDGEC3.75CM) 


Ci 


CELLS 2125 -X% sizs Geounp 
PILE FOR PLOQIENG USING PLOnEGur EXEC 
MOO 55 obpeo Sal ptno. Sadly oleal 


Ved 0050, 0, 05-1 25,000 -100 
200,6, SS Oe oe eo eo 
2,550,0,0,0, 225,05 200ecor 
5,5,0,0,0,—.12550,0,200.201 
S004, ..5, =: 5 ly. oe O noe ook 
457,060,050, 0,0s-.125,500.5035 
S003). 55 —20/ Spl Ue 57 Say ol 
$,2,0,0,0,0,.125,0,4600.906 

SOG sG5 Sr SpSo Beds - S/S, 0,31 - (o} at 
Gen, 0,0, -.lc5, 0707 S00, 505 
0,110 

996 ye > 07 O pee So Dipl Oe 
7,3,0,0,90,0,0,0,099.099 
256505 =. 5e15 Ope ool 

Snee oSinWse lama Sn Un la o Wal 
lvo.0.005,0,1,0.005,071.6, ,016 
0535,0,0,10,0,0,0, 001 00) 
O70,0,0;,0,—.5/5-0,0,-001 002 


Data File Geometry of 6cm Monopole at Edge 
of Wire Grid veo 
191 


FILE: NEOGIt DATA Bl 


CM WIRE GRID 10 CM. LONG WITH MONOPOLE 6 CM AT CORNERC5.3CM ON DIAGONAL) 
Cit "CELLS 2125 % ples Geouns 

CE FILE FOR PLOTTING USING PLOTUG RE Ve tee 

G} L00;6, 5, = 05,1505 =.5,707 201 


Git L785 050,05 Oy =o Ue 100 

GW 200,8,.5,- ils obs —oeln Uy Ol 

GHZ, 5, 8,05 0705. leo. 0s 2uUezol 

Gi4 Sy, OO Oh cbr Or Oc Z01 

Gh S006 S> .Sya eo le el Lee 01 

Gi4 G,7,>0,0,050505> 7lco, 500 .s05 

GY GUC) S,..5,-257 5.0, 0, -~57 So, 01 
iM 5,2,0,0,0,0,.125,0,400.4904 

GH C0 iS oS enn en llA cei! Gin Ue tbh oul 
GM 6,2,0, 0,0), — 12570 Ol 5 00 S105 

GX 0,110 

GW 99,65 .50,50,05.5,0,1,..01 

Gi 7,3,9,0,90,0,0,0,099.099 

GN Pelee Wirwadeno prs bere Vics eosin ade a 

GW SyO pss Ue = oo Us On 

GH 165.0, 005, 05150 .005,0,1.65.006 

Gi4 0,35,0,0,10,0,0,0,001.001 

GM) 0505.0,05 0,257 53.57 570,001, 0G 

GE 

EN 


3. Data File Geometry of 6cm Monopole at Corner 
of Wire Grid Box. 
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PART B: 


NEC INPUT DATA FILES FOR NEAR ELECTRIC FIELDS. 


File: Wl DATA Bl 
CH vilRE Grldetg sc, LONG 
Cit Celis Wieser 625: GROUND 
Ce 
Gil LOO 54; See ole ops CL 
Gi4 Lara eee ace, LOO. 1 00 
Gil 20032 Spee oe oo OL 
Git oe ss OO De Une cop OeeUOe eu) 
Gi4 Sia Cag Choy 1 te Greeny 0a eal a| a ea 
Gil SOO hey oSpHaSnls oSo Wp lo a ul 
Gi Sr AWaG sly Ws bl, —, LAE 5 Solo} , Sols 
GH GO0GG. coe Oe tr Os Sool oo Perel 
Gi4 5,2,0,0,0,0,.125,0,400.404 
GH S00 -S7 457 Ss —1 Oe ores 2b 
Git 6.250.005 bcs. Ure, eo. oS 
A 0,710 
GSO 70g. 
Gel 
PR 070,020, 2000 
GH 1 
WG 
NX 
Cil 6 Cy MONGPOLE Al CENTER ®GF NIRE GRID FREQ=1 GHZ 
= NEark E=rPPebp on A-axkisS. CCONTOURS 30 POINTS) 
GF 
Ghl 2716 05 = 51, Ol, Sp ehn o Ol 
Gi. Specie ei nla eS Sra dhs 3 191 
GH 1S Sin [0 ale A I) ’ sn o UIE 
Gs 0,835.1 
Ge 
Ba UG lol Clore 
Pare oe 
NE 0, 350,15.35050. 0166,0,0-115,.0.),.0, 0.1 
XQ Q 
EN 
1. Data File 6cm Monopole at Center of Wire 
Grid Bex (fer Contour Plot of Magnitude). 
Po IE SSP Bi DATA Bl 
CM WIRE GaiD 10 cm. Lone 
so CEULS lZs 5%, el 25 Cr OUrnD 
GW 10059) 55-3557 0 pe eo 
Gil 1,7,0,0;, 0,0 70> —-.2s.00e oe 
Gi 200733. 5 = 57 le eee ee 
Gil 2,35,0,0,0,0,. 125,50, 200,201 
Gi4 55520 70,0, > leo, 0,U aco col 
GH 50079, 553-5 olor OL 
GM 4$,7,0,0,0,0;0,—.125,500.505 
GH 4600,4,%5/3257 551 0, -1575,1, Ol 
Git $,2,0;50;0,0,.125, 0,400. 405 
CGAC500,5, .57S7- eo le or Oe o Lt 
Gil 6,2,0,0,0,—=.125,0,0,500 505 
Aor L110 
SSO, 1 
GE. 
FR 0,0,0,0,1000 
Git) 
viG 
NX 
C1 6 CM MONOPOLE AT CENTER OF WIRE GRID FREQ=1 GHZ 
CE NEAR E-FIELD GON X-AKIS “PHASE 100 POINTS) 
GF 
GW 236; 05 — 255 15 05 oe eo 
Gh S. 00 «5) 0s ls oe SoU pelo 
GH VW, S,0; 05). 0.0. .05000 bo 
G5 0,0,.1 
GE 
~ Op lel, Cyl .u 
PL 2,2,4%,1 
HE 0,100,.1 ,100;,020166,0, 0.115) 05029,0.0.929 
XQ 0 
EN 
2. Data File 6cm Monopole at Center of Wire 


WIRE GRID MODELING 


Grid (lor Conrour seleteot Phase ). 
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FILE: Wess DATA Bl 


HIRE GRID 10) Che Jeune 

CELLS .125 K% .l2sveecunT 

LUG; 4.5 = slo ea By 30)! 
ly? 70>0,0,0,05 =e co, Core ue 

1 200-8. 8,2 (504 Mab mees Jomeclt 
2r5>0;70,070,.225,07 20Uu coe 
3,3,0,0;0,=.125,0,0,e007cu0) 
500,46, .5,=<5)15 8205. ue 
4,7,0,0,0,0;, 0, 772257 Gurus 
400,46, .5,~.35755 ),0,-—-57 5, Ol 
§,2,0,0,0,0,.125,0,4900. 40% 
500,55 .3/755~.55ls oo Math » Wal 
6,2,0,0,0;—-.12570, 0,200,505 
0,110 

0,0, .1 

il 

0,0, 0, 0,1000 

a 


6 CM MONOPOLE AT EDGEC3.75 CM FROM CENTER) OF WIRE GRID FREG=ieeae 
NEAR E-riBlED ON X-AXiS “CONTOURS 20) POINTS) 


PERS Ne a Sip db A Cle Syren a 

Spo c Big Ogee I OF eae 

1757 257 5 Uso ous ae 6,.016 

0,0, 1 

Opel ts 0 

Zp er ore 

0,750.1; 30, -002095 02.07 1) Ss Os Ole 
0 


Data File 6cm Monopole at Edge of Wire Grid Box 
(tor Contour Plotwot Magnitude), 


FILE: WEDPH! DATA Bl 


cM 
Cm 
CE 
GN 
GM 
GN 
Gi4 
Gi 
Gi 
GM 
GW 


WIKRE GRID 10 CM. Cone 
CELLS 2.125 X% 2i25 Ground 


100; G,. 5 55 1s 01 ed 5, Nee 
Lind. 0. (noes lZo> TOOe MA 
ZU Se sa0 Sm oie tae, anos SU OL 
ran a Ah s(0l Aiilp (O)Be 125,0,200.201 
3757070, 05-. T25s00 0, 200.201 
S005 9). 55>. Syl 5; Orr eee): 
4,7,0,0,0,0,0,-, 1257500505 
G00,4>5 .95 —. 375s 150,257 5. Agee 01 
5,2,0,0,0,0,.125,0,400.404 
800545337 55->. 55 lees oo ea) 
6,2,0,0,0,~.125,0,0, 500.505 
0,110 

05 0s%.2 

Ag 

a »0,0,0,1000 


6 CM MONOPOLE AT EDGEC3.75 CM FROM CENTER) OF WIRE GRID FREQ=]eGHe 
NEAR E-FIELD ON X=AXIS “GPHASE 200 sFOLNTs2 


“ 
a 


00,-0.0209,0,0.115,9 .02950, 02022 


Data File 6cm Monopole at Edge of Wire Grid Box 
(for Contour Plet sor eehasces- 
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Fab e: WCCO] DATA Bl 


CM 


WIRE GRID 10 CM. CONG 


ae GCELES .125 xX .12536R0Cno 
c 
GW 0074s «5, —~.5>1,05>. See 
Gil, ?,0,0,0,0,0,-.125, 100.100 
Gil 20075) 59-1. Ox lp wopeoe Uo 
Garcs>s,0,0,0,0,.125,0,200. 201 
Grips, 2,7 0,0,0,-.125,0,0, 200-2071 
GW 500; 4, sap =e alla o Sele pe 01 
Me9,7,0,0,0,0,0,-.125,500).3035 
GW SOG 4s. Os 7257 So lp Ue —nar els 01 
GM 5,2,0,0,0,0,.125,0,400.404 
en 500,49, .375,-.5,1,.575.0,8, 01 
t G,270.0,0,—-sl25,0, 00,5005 505 
m™ 0,110 
ea, 0, .1 
GE l 
FR 0,0,0,0,1000 
GN l 
WG 
CM 6 CM MONOPOLE AT CORNERC5.3CM ON DIAGONAL) OF WIRE GRID FREQ=1 GHZ 
CE NEAR E-FIELD ON X-AXIS (CONTOURS 30 POINTS) 
GF 
Bice e, 0,—.521,0,.5,1, 201 
GH Oop Bhp OF e— Ge OR ’ ol 
Baglo 5.5755 .57 5, 1a. 3750. $75, 1.365.016 
Gs 0,0;.1 
GE 
reo, 1,>1,0,1,0 
meee, Dl 
HEotns 0,2, 50,0.0541,0.0375, 0-115 021507002 
XQ 0 
EN 
5. Data File 6cm Monopole at Corner of Wire Grid Box 
(for Contour Plot of Magnitude). 
Pole: WCOPHIL DATA Bl 
CM WIRE GRID 10 CM. LONG 
ot CEEES «125 X .i25 GROURD 
CHGS, 55 7. Sele Os edo ks «0 
Giel,7,0,0,0,0,0;,—-.125,100. 100 
GN 20085 .53 7. Dp less 7 Dp 0s Ue 
Giee, os, 0,0,0,0,.125,0,;200.201 
GM Sp So Up WnW neon Wnts Zu) 62a 
GW SHH 545 -Spes ba ly Pas Wry gran lier ca iI 
Gi4 Sea 0, 0, 0, 0,>. eo, 00.00 
GW Oa 255-57 Oy Le pm ot Oe ky 01 
GM 5,2,0,0,0,0,.125,0,400.404 
Cimo0e 4, 91 5)7.5+1,.57550), 1,01 
GM Gc, 0; O50 e— 1255050; 500.505 
GX 0,110 
G5, 0,0, .1 
GE l 
FR 0,0,0,0,1000 
GN 1 
i 
N ‘ 
CM 6 CM MONOPOLE AT CORNERC5.3CM ON DIAGONAL) OF WIRE GRID FREQ=1 GHZ 
ae NEAR E-FIELD ON X-AXIS CPHASE 100 POINTS) 
ONG 270, 07-.5,1,0,.5,1,.01 
Giibsgeg-5,0,1,=.5,0,15.01 
Gel, 5,.3575,. 375, oe 5154. 51 Op eo OLS 
GS Geet 
GE 
EX id AU a 
PL 2,2,4,1 
Meno, 20071,100,0.05461,0.0375,0.115,0.029,0,0,029 
an 
0 
EN 
6. Data File 6cm Monopole at Corner of Wire Grid Box 


(for Contour Plot of Phase). 
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